Changing universe model with apphcations 
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Although the current galaxy models yield calculations consistent with much of the data, many 
irregularities exist, exceptions have been found to the current models, the ACDM model apparently 
fails on galaxy scales, dark matter remains elusive, the phenomena at the center of galaxies are only 
beginning to be addressed, and many observations and empirical relationships are unexplained. A 
changing universe model (CUM) is proposed that posits the stuff of our universe is continually 
erupting into our universe from sources at the center of galaxies. In this first test of the CUM, a 
single equation that describes the rotation velocity curve of spiral galaxies is derived. The equation is 
also used to correlate the measured mass of the theorized, central supermassive black hole with other 
galaxy parameters. The equation adds one parameter for each mature galaxy and one parameter for 
each particle (atom nucleus and smaller) species to Newtonian dynamics. Rising, flat, and declining 
rotation velocity curves are explained without unknown dark matter. Ten previously unknown 
parametric relationships are discovered. Also, the CUM integrates several heretofore unrelated 
observations. The CUM model is in an early development stage and has already predicted new 
relationships among galaxy parameters. 

PACS numbers: 98.80.Bp,98.62.Ai,98.62.Dm, 98.62. Js 



I. INTRODUCTION 



The discrepancy between the Newtonian estimated 
mass in galaxies and galaxy clusters and the observation 
of star and gas kinematics is well established. Currently, 
the possible explanations for this discrepancy are that a 
large mass Mdm of dark matter (DM) and dark energy 
exists (the ACDM model) or that Newtonian gravity is 
modified on galactic scales. The current ACDM model of 
the structure of the universe appears successful on large 
scales pl , but galactic scale predictions disagree with ob- 
servations 0, 0, 13 ■ The most popular modified gravity 
model is proposed by the Modified Newtonian Dynam- 
ics (MOND) model ^ and references therein]. MOND 
suggests gravitational acceleration changes with galactic 
distance scales. MOND appears limited to only the disk 
region of spiral galaxies and appears to have two falsi- 
fiers m. How MOND may be applied to cosmological 
scales is unclear. However, MOND may represent an ef- 
fective force law arising from a broader force law 0| . The 
evidence suggest the problem of a single model explaining 
both galactic scale and cosmological scale observations is 
fundamental 0. 

An alternate explanation is the possibility there ex- 
ists another force originating from the center of galax- 
ies. This proposal, called the changing universe model 
(CUM), posits the stuff of our universe is matter and 
a matterless scalar p whose gradient exerts a force Fs 
on the cross section of matter. The Fg is outward from 
galactic center. Thus, Fg is a repulsive force in a galaxy 
relative to gravity. However, it differs from a gravity type 
force since it acts on the cross section of matter and it is 
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matterless. 

The rotation velocity v (kms""'^) of a particle in the 
plane of a spiral galaxy (hereinafter galaxy) reflects the 
forces acting on the particle. An explanation of v'^ as 
a function of radius r (kpc) from the center of a galaxy 
(RC) requires a knowledge of galaxy dynamics and an ex- 
planation of radically differing slopes. The v and, thus, 
the RC is measured along the major axis. Although v is 
non-relativistic, the result of calculations using RC mod- 
els are used in cosmological calculations. Thus, the RC 
is an excellent tool to evaluate galaxy models. Further, a 
complete explanation of the RC must also explain black 
holes and observations of the center of galaxies. 

Battaner and Florido and Sofue and Rubin Q pro- 
vides an overview of the current state of knowledge of 
RC's in the outer bulge and disk regions of galaxies. The 
particles most often measured in the disk region of a 
galaxy are hydrogen gas by Hi observation and stars by 
observing the Hq, line. The particles being measured in 
the inner bulge region are stars by observation of the of 
Hq, CO, and other spectral lines. Also, the RC differs 
for different particles. For example, although the Hi and 
Hq RCs for NGC4321 H in Fig. □ differ in the outer 
bulge, they approach each other in the outer disk region. 
Most models of RC monitor the Hi in the outer bulge 
and disk regions and ignore the differing Hq RC. Also, 
when drawing the Hi RC, a straight line is often drawn 
through the center. This practice is inconsistent with 
stellar observations ^3]. 

The observation of rising RCs in the outer bulge and 
rising, flat, or declining RC's in the disk region of galax- 
ies poses a perplexing problem. If most of the matter 
of a galaxy is in the bulge region, classical Newtonian 
mechanics predicts a Keplerian declining RC in the disk. 

Currently, a convention on how to calculate error (e.g. 
random or systematic) in RCs is lacking in the literature. 
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FIG. 1: Plot of the Hi (filled diamonds) and Ha (filled 
squares) RC for NGC4321. 



One popular method is based on the assumption of sym- 
metry. This practice of averaging the v of the approach- 
ing and receding sides of the RC is ignoring important 
information [Tll |. Asymmetry appears to be the norm 
rather than the exception [ij. The RC of all nearby 
galaxies for which the kinematics have been studied have 
asymmetry ^3 ■ larger distance, objects appear to be- 
come more symmetric as they become less well-resolved 
(Hi resolution is typically 15 arcsec to 50 arcsec) [ll|. 

Hodge and Castelaz jljl found a correlation between 
the size and distance of neighboring galaxies and the 
asymmetry in the disk of the Hi RC. The galaxies studied 
with minimal asymmetry also had rising RCs. Therefore, 
rising RCs are intrinsic to galaxies. 

A strong correlation between the core radius Rc and 
the stellar exponential scale length Rd was found by Do- 
nato and Salucci ^3 ■ The Rd is derived from photomet- 
ric measurements. The Rc is obtained from kinematic 
measurements. If this relationship holds for both high 
surface brightness (HSB) galaxies, low surface brightness 
(LSB) galaxies, galaxies with rising RCs, and galaxies 
with declining RCs, modeling such a relationship is very 
difficult using current galaxy models. 

Chez et al. E3 Ferrarese and Merritt 0| have 
observed Keplerian motion to within one part in 100 in 
elliptical orbits of stars that are from less than a pc to a 
few 1000 pc from the center of galaxies. 

The orbits of stars within nine light hours of the 
Galaxy center indicates t he p resence of a large amount 
of mass within the orbits "11111111111 El- To achieve 
the velocities of 1300 kms'^ to 9000 kms'^ |l3 and 
high accelerations 0, there must be a huge amount 
of very dense particles such as millions of black holes, 
dense quark stars j20l and references therein], and ion- 
ized iron |2lj | distributed inside the innermost orbit of lu- 
minous matter. The mass M, (in Mq) within a few light 



hours of the center of galaxies varies from 10^ Mq to 10^^ 
Mq |23| . The M, can be distributed over the volume 
with a density of at least lO^^ Mq pc'^ El El El 11 . 
The orbits of stars closest to the center are approximately 
1,169 times the Schwartschild radius of the theo- 
rized supermassive black hole (SBH) E3- However, such 
large mass crowded into a ball with a radius of less than 
60 AU must either quickly dissipate, must quickly col- 
lapse into a SBH 25, 26], or there must exist a force in 
addition to the centrifugal force holding the mass from 
collapse. Mouawad et al. |23] suggested there is some 
extended mass around Sgr A. Also, models of supermas- 
sive fermion balls wherein the gravitational pressure is 
balanced by the degeneracy pressure of the fermions due 
to the Pauli exclusion principle are not excluded ^2S\. A 
strong repulsive force at the center of galaxies would ex- 
plain the shells of shocked gas around galactic nuclei p9l 
page 59 5] |3C1|. the apparent inactivity of the central ob- 
ject 1^ I32ll33l 134] , and the surprising accuracy of rever- 
beration mapping as a technique to determine M, [35l] . 
However, discrepancies have been noted among the meth- 
ods used to measure M, [sl [sl and references therein] . 

The first published calculation of the slope of the M, 
to velocity dispersion tJv (in km s~^) curve (M, — ctv re- 
lation) varied between 5.27±0.40 [33 and 3.75±0.3 [13 . 
Reference [H and references therein] suggested the range 
of slopes are caused by systematic differences in the ve- 
locity dispersions used by different groups. However, the 
origin of these differences remains unclear. 

Ferrarese |39j found the ratio of the M, to Mdm of 
the theorized dark matter halo around galaxies was a 
positive value that decreased with halo mass. However, 
if the intrinsic rotation curve is rising the effect of 
the force of A/dm in the equations implies the effect of the 
center object must be repulsive. Such a repulsive force 
was called a "wind" by Shu et al. 40], Silk and Rees j^. 
The "wind" (a gas) exerted a repulsive force acting on 
the cross sectional area of particles. Therefore, denser 
particles such as black holes move inward relative to less 
dense particles. 

A multitude of X-ray point sources, highly ionized iron, 
and radio flares without accompanying large variation at 
longer wavelengths have been reported near the center of 
the Milky Way [H il IH IH iHl HI HI . 

Mclure and Dunlop ^3 found that the mass in the cen- 
tral region is 0.0012 of the mass of the bulge. Ferrarese 
and Merritt 116| reported that approximately 0.1% of a 
galaxy's luminous mass is at the center of galaxies and 
that the density of SBH's in the universe agrees with the 
density inferred from observation of quasars. Merritt and 
Ferrarese |45| found similar results in their study of the 
M, — Ctv relation. Ferrarese |3^] found a tight relation 
between rotation velocity Vc in the outer disk region and 
bulge velocity dispersion CTc (vc — CTc) which is strongly 
supporting a relationship of a center force with total grav- 
itational mass of a galaxy. Wandel ^3 showed the M, 
of AGN galaxies and their bulge luminosity follow the 
same relationships as their ordinary (inactive) galaxies. 
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with the exception of narrow Une AGN. Graham et al. 
I47L liM li^ found correlations between M, and struc- 
tural parameters of elliptical galaxies and bulges. Either 
the dynamics of many galaxies are producing the same 
increase of mass in the center at the same rate or a feed- 
back controlled mechanism exists to evaporate the mass 
increase that changes as the rate of inflow changes as sug- 
gested by Merritt and Ferrarese 35]. The RCs imply the 
dynamics of galaxies differ so the former explanation is 
unlikely. 

In accordance with the Principle of Fundamental Prin- 
ciples (see Appendix 0, a search was made for a physi- 
cal process in the Newtonian physical domain that may 
model and explain the observed data of spiral galaxies. 
Such a system was found in the fractional distillation 
process (FDP). A FDP has a Source of heat (Space en- 
ergy) and material input (matter) at the bottom (center 
of galaxies) of a container. The heat of the FDP is identi- 
fied with a constituent of our universe that is not matter 
and is identified as Space (in units of "que" , abbreviated 
as "q" , - a new unit of measure) whose density p is a po- 
tential energy. The term "Space" will mean the amount 
of Space in ques. In the FDP, temperature decreases with 
height from the Source. The opposing forces of gravity 
and the upward fiow of gas acting on the cross section of 
molecules in an FDP cause the compounds to be lighter 
with height and lower temperature. Heavy molecules 
that coalesce and fall to the bottom are re-evaporated 
and cracked. The heat equation with slight modification 
rather than gas dynamics modeled the flow of heat (Space 
energy). Temperature decline of the FDP was modeled 
as dissipation of matter and Space energy into three di- 
mensions. Thus, the Changing Universe Model (CUM) 
was born. 

This paper posits that there is a Source of strength e 
at the center of each galaxy. The Source is a source of all 
matter and of Space. The p forms a scalar potential field. 
Energy is conserved in the emission of matter/energy and 
Space from the Source. The unit of measure of e is in the 
amount of que released per second which is proportional 
to the amount of mass released per second, Cm in units of 
M0 s"^ The apphcation of the CUM to the RC of galax- 
ies explains observations and empirical relations. Also, 
ten new relationships among galaxy parameters were pre- 
dicted and confirmed. This, the first test of the CUM, 
successfully explains RC observations and suggests an al- 
ternative to DM. The CUM, at this early stage of model 
development, is bringing together heretofore apparently 
unrelated observations. 

The object of this article is to develop and test the 
CUM by suggesting a galaxy model consistent with the 
CUM, by explaining observations of the RC and the cen- 
ters of galaxies, and by deriving relations among galaxy 
parameters. In section ^ the CUM model is developed 
from the fundamental principles, which are stated in Ap- 
pendix ^ An equation for the RC is developed from 
the energy equations of the CUM, applied to galaxies, 
and compared with galaxy observations in Section IIIII 



The posited description of particles is used to explain 
several phenomena involving observations of the center 
of galaxies. Other parameter relationships are presented 
in Section Hvl Section Ivl shows the interger relatinships 
developed in the previous sections are not random. The 
results are discussed in Section IVIl Section IVlII lists the 
conclusions. 



II. CUM MODEL 

The assumed Fundamental Principles are stated in Ap- 
pendix^ Throughout the text, capitalizing the first let- 
ters of the name of the Fundamental Principle denotes 
these principles. 

The Principles of Change and Repetition requires that 
the three dimensions (3D) be created from two dimen- 
sions (2D). The creation of 3D from zero dimensions (OD) 
is not only unlikely but forbidden. The universe begins 
with a Change that may be a rupture of OD. The Change 
itself is a beginning of time. The OD ruptures to cre- 
ate one dimensional (ID) points and Space between the 
points on a line. The OD is a Source of ID and con- 
tinually ejects Space into one end of the ID line. The 
other end must expand into the void. Thus, the ID line 
is continually growing. The ID line has only one direc- 
tion. Therefore, the ID line is straight. The ID energy is 
Tit — pil where ri is the point tension, t is time elapsed 
between point eruption, pi is the line Space density, and I 
is the length between points. The first transition creates 
length and time 

I = ct, (1) 

where c = Ti/pi is the constant of proportionality. 

When points in ID rupture, the 2D Space is inserted 
into 2D at many points along the ID line. The rupture 
can have no component remaining in ID. Thus, a right 
angle is defined and the metric coefficients equal one (a 
Cartesian plane). Accordingly, 2D is flat. The 2D en- 
ergy for each transition is T2I — P2S where T2 is the line 
tension, / is the line length, p2 is the area Space density, 
and S is the area. 

A 2D Sink is at a place in 2D that has Space tightly 
contained in a line. The line is a one-dimensional border 
in 2D. If a line with a "line tension" compressing Space is 
deformed perpendicular to the line, an arc of height h is 
formed. The geometric constant of the line is its length 
I. If a condition oi I ^ wh is reached, a bubble in 2D is 
formed. The end points of the line are coincident. This 
rupture causes Space to be removed from the 2D plane. 
This 2D Sink becomes a Source in 3D Space. The line 
becomes an object in 3D called a hod. The hod is the 
smallest piece of matter from which all other matter is 
constructed. The hod is supple so that its shape may 
change in 3D. The edge of the 2D plane expands into 
the void. Hence, the distance between the points that 
rupture into 3D is increasing. The expansion is recur- 
rent and is in 2D, not 3D. Because the rupture into 3D 
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is orthogonal to 2D, 3D is flat (Cartesian three space di- 
mensions), also. The I becomes area a. The Space area 
in 2D becomes the volume V of Space density p in 3D. 
One que is the amount of Space released by one hod in 
the transition. 

Since nothing exists in 3D except because of the tran- 
sitions, the energy for a transition is 

™ = (2) 

where pi and Vi are the initial p and V at the transition, 
respectively, and r is the proportionality constant called 
surface tension. 

Similarly, corollary II of the Fundamental Principle im- 
plies a 3D Sink may be formed. 

Succeeding ruptures at the same point adds Space and 
hods to the universe. Since we observe there is Space and 
p in our neighborhood, the Space and hods then flow from 
the Source and creates the galaxy around the Source. 

The total energy transferred into 3D from Eq. is 

NtTa = NtpiVi, (3) 

where A't is the total number of hods in the universe. 

In the observed universe, the 2D plane is common for 
all Sources. As the Space from a Source expands into the 
void, it encounters other Space from other Sources. Since 
all Space originates from one OD Source, the merged 
Space is causally correlated and coherent. However, 2D 
Space is distinct from 3D Space. 

In 2D, the surface was a tightly bound line on the 
outside of Space. In 3D, Space is on the outside of a 
surface. The r on the 2D Space was high enough to 
contain the Space. Now r acts to attract 3D Space. Since 
the hod formed in the extreme pressure of the Source, the 
T maximum ability to act on Space exceeds the pressure 
of Space on the hod away from the Source. Since r and a 
are constant, as the extent of the volume V increases so p 
must decrease. Therefore, the hod applies its energy, ra, 
on the Space immediately surrounding the hod. Since 
the hod contained the Space before release and Space 
was at a much higher p than is surrounding the hod in 
3D, all the hod's tensional energy is exerted on Space. A 
hod is a 2D discontinuity in the continuity of 3D Space. 
These discontinuities are surfaces in Space with an edge 
not containing loops. 

Hods are discrete. The mathematics of integer num- 
bers can then be used to count hods. In addition, hods 
can combine, one with another. The number of hods and 
the structure of the combinations of hods create the dif- 
ferent particles and forces, except gravity, of the universe. 
Therefore, the language of hods is integer numbers, struc- 
tures, counting, subtraction, and summation. 

Space is continuous, perfectly elastic, and compress- 
ible. The Space continuously fills the region it occupies 
except for the hods. Transformation of Space over time 
is continuous. Hods and the outer edge of Space as it 
expands into the void are discontinuities in Space. Space 
really can meet the mathematical criteria of considering 



an arbitrarily small volume. Therefore, the language of 
Space is real numbers, tensors, field theories, differentia- 
tion, and integrals. 

The three interactions of concern are that each hod 
exerts a force on Space, Space exerts a force on the hod, 
and Space exerts a force on adjacent Space. A hod is 
not attracted or repelled by other hods through a field 
or other action at a distance device. Since r can act only 
perpendicular to the hod's surface, the force Space exerts 
on the hod is perpendicular to the surface of the hod. 

There are three dimensions in our universe. Counting 
is a characteristic of hods. Therefore, that our universe 
is three-dimensional is a characteristic of the hods and 
not of Space. 

There are two concepts of distance. Distance can be 
by reference to a fixed grid. Traditional methods have 
used Space as the fixed grid. In this model, distance is 
a multiple of the average diameter a hod. Since a is a 
constant and the hod may be deformed, the hod average 
diameter will be considered a constant throughout the 
universe beyond the Source and in moving and acceler- 
ating reference frames. The measure of distance will be 
a "ho" , which is defined as the average diameter of the 
hod. The CUM uses ho distance to form a grid in a fiat 
universe with 

lho = 2(-)i/2. (4) 

TT 

Another concept of distance is the amount of Space 
between two surfaces. By analogy, this type of distance in 
gas is the amount of gas between two surfaces. Increasing 
the amount of gas between the two surfaces increases the 
distance. One way to do this is by movement of the 
surfaces. Another way is to add more gas from a source. 
If the medium is compressible, the distance measured 
by the amount of gas varies with the change in density 
of the medium. The CUM model allows Space to be 
compressible. Therefore, the amount of Space between 
two surfaces can change if the p (que ho^'^) of the Space 
between two surfaces changes. 

The terms distance, length, area, and volume will refer 
to ho measurements. The unit of measure of distance is 
usually expressed in km, pc, kpc, or Mpc. If the mea- 
surement is taken by means of angular observation or 
time-between-events observations such as using Cepheid 
variables that do not depend on the amount of Space, 
the units of pc are proportional to the ho unit. Since the 
proportionality constant is unknown, the units of pc will 
be used. However, if the measurement such as redshift 
depends on the amount of Space, the units of pc should 
be used with caution. 

Since the hods are discrete, e{t) is a digital function in 
time. The hod exerts a tremendous r to hold the Space in 
its 2D surroundings. Upon transition, the Space immedi- 
ately around the hod experiences an instantaneous addi- 
tion to the amount of Space in the 3D universe. This step 
function creates a Shockwave in the surrounding Space. 
Space then flows out into the void causing the V of Space 
to increase and a p variation in the V. 
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If Space has no relaxation time and no maximum gradi- 
ent, the eruption of Space from a Source would spread out 
instantly so that p would be virtually zero everywhere. 
Since we observe Space in our region of the universe, the 
Anthropic Principle implies Space flow must have a re- 
laxation time and a maximum p gradient ( Vm/f) that is a 
constant of the universe. Bold face type denotes a vector. 

Since Space has a relaxation time, the mirror property 
in hods is inertia i per a. A hod is an object with an 
area a, a surface tension per unit area r, and inertia per 
unit area t. The a, the magnitude of r, and the l are the 
same for all hods and are universal constants. 

The T and t furnish a resistance to deformation on the 
hod. Since a is a constant, a deformation of a hod causes 
the edge to change shape. 

In the region around an active Source, Vp = Vmp. 
Space expands into the void where p — and Vp = 0. 
Therefore, at some radius Vgc from the Source, p becomes 
another function of time and distance from the Source. 
At Tsc, VjO must be continuous, Vp = V,„p and V^p > 
0. If e{t) changes, then changes. 

As the Space increases V , Space increases between the 
hods. The hods move and attain a velocity ddij /At where 
dij is the vector distance from the i"^ hod to the j"^ hod. 
The hod's l inhibits the movement. Thus, a feedback 
condition is established wherein the ddij/dt is balanced 
by t. Therefore, the hods movement is proportional to 
the total inertia La of the hod and to the velocity of the 
hod. The momentum pij between the of the i*^ and j^^ 
hods is defined as 



Pij 



Kp La 



~df 



(5) 



where Kp is the proportionality constant and = means 
"is defined as" . 

The energy Ty expended by the Space expansion of V 
upon the hod is 



Tij = / Pij •dd = KpLa / — 



/- 


ddij 


/ dt 


dt 



dt. 



(6) 



between hods, the energy the hods exert on each other 
must increase. Since the p term of Eq. ^ does not apply 
to the hod, the loss of energy must be due to ra. Call 
this energy between the of the z*** and j^^ hods the poten- 
tial energy Uij. By the Principle of Negative Feedback, 
Uij is proportional to ra and is a function f{dij) of dij 
between hods 



Uij = ~Ku ra J {dij ) , j 



(9) 



where is the proportionality constant and the nega- 
tive sign is because the hods exert an attraction to each 
other. 

The total Uij energy Uj exerted on the j*'' hod by all 
other hods is 



Uj= J2 



(10) 



The j hod's ta inhibits the changing pij . The inertia 
of the hods results in work W done on the hod in time 
dt 



Wj = - 



r d ^' 



d^i 



Wi = -K^ La ( 4 " ^^^^ - db- 



dt^ 



(11) 



(12) 



where subscript "b" and subscript "e" indicate the be- 
ginning and ending positions in the time period, respec- 
tively. 

Therefore, the total energy equation in 3D is Eq. Q 



with Tj and Uj is 



Nt 



Thus, 



N,Ta+ +Uj) ^ /// pdV = Q. (13) 

i^jj=l JJJSpaco 



Tij = Kt Kp La 



ddi 



dt 



(7) 



where Kt is the proportionality constant. 

The Tij is a scalar quantity and is the kinetic effect 
of the i*'' hod on the ^od. The total energy Tj 
exerted on the j**^ hod by all other hods is 



T 



Nit) 

3 — X/ ^'J' 

i^j i=l 



n. 



(8) 



Because the amount of energy added by the movement 
of La must be the negative of energy due to position and 
by the inertia which inhibits hod movement, the Prin- 
ciple of Negative Feedback applies. As Space is added 



Implied in Eqs. © and is 

dd. 



Nt Nt 



dt 



= E E^"™/('^^^)- 

(14) 

"Flatness" is a fundamental condition of the creation 
of matter and Space. Therefore, the Cartesian coordinate 
system with four equal angles to define a complete circle 
is the form of the "flat" measuring system. 

Note, dij is between hods. The process of summing 
over all hods provides a reference frame similar to Mach's 
principle. Also, since there is distance between galaxies, 
the Uij and Tij of hods in different galaxies includes a 
heritage from ID and 2D. 
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A. Energy continuity equation 

Consider an arbitrary volume AV^ bounded by a sur- 
face As. Space will be in AV. Hods, Sources, and Sinks 
may be in AV. Equation {Tj) implies = Tji. Equa- 
tion lO implies Uij — Uji. Therefore, the total energy 
in AV is 

= Nra + / / / pdV 

J J J AV 

N 

+ Y.^T, + U,+W,) 

fc=i 1=1 



(15) 



where A'' is the number of hods in AV, the integration 
is over AV, is the number of Sources in Ay, iV^ is 
the number of Sinks in AV, ek is the strength of the fc**^ 
Source, and rji is a negative number which is the strength 
of the Sink. 

The classical development of an energy continuity 
equation includes the assumption that energy is neither 
created nor destroyed. In the CUM, an arbitrary vol- 
ume may include creation into our universe of energy 
(Sources) or the destruction from our universe of energy 
(Sinks). 

Following the usual procedure for the development of 
an energy continuity equation such as the heat equation 
for the temporal evolution of temperature; remembering 
hods have no volume and, therefore, the region can be 
made regular and simply connected by mathematically 
drawing cuts around the hod discontinuities; and using 
the Principle of Minimum Potential Energy and the Prin- 
ciple of Feedback and Minimum Action gives the integral 
form of the energy continuity equation 



C—E^dV 
at 



S 1 1 Vp» nds 

N 



St V 



N 



N 

E 



' nds 



• nds 



[TaNtuw{t)vn] •nds 



N„ 



(E^fc + E^')d^' (16) 



and is assumed to be independent of position or direc- 
tion. S' is a factor that slows the flow of Space. The idea 
of their being some kind of viscosity or friction associ- 
ated with Space would require the expenditure of energy. 
Since Equation H13() lacks such a term. Space energy dis- 
sipates without loss of total energy. Therefore, there is a 
relaxation time associated with Space. This means that 
the stress created by the sudden entry of Space upon 
a transition must relax asymptotically with time. The 
compression of the transition is dissipated by the flow of 
Space into the void. (2) The Su is a constant of propor- 
tionality. Call Su the potential conductivity. S^i is the 
ability of energy to flow and is assumed to be indepen- 
dent of position or direction. (3) The St is a constant 
of proportionality. Call St the kinetic conductivity, 
is the ability of kinetic energy to flow and is assumed to 
be independent of position or direction. (4) The S'w is 
a constant of proportionality. Call 5^ the work conduc- 
tivity. S'w is the ability of work energy to flow and is 
assumed to be independent of position or direction. (5) 



The iVt„w = N 



EUiT, + U, + W,)\t,,u„w, is the 



number of hods which leave the AV with the energy as- 
sociated with each hod and \Tj,Uj,Wj means that Tj , Uj, 
and Wj are constant with respect to the differentiation. 
(6) The Vn is the velocity vector of the A^tuw 

Applying Gauss's divergence theorem to the surface 
integrals of Equation and noting that the volume 

could be chosen to be a region throughout which the 
integrand has a constant sign yields the differential form 
of the energy continuity equation 



d£:, 



N 



E 



c 



N 



E 

+ V • [TaNtnwit)Vn] 

+ C-i(^e, + ^,7/), 



c 



N 



E 



(17) 



fc=i 



1=1 



where = S/C and is called the diffusivity constant of 
Space. 

The presence of Sources and Sinks leads to a non- 
homogeneous equation. Note the similarity of this equa- 
tion to the diffusion equation for the temporal evolution 
of the concentration of chemical species in a fluid and the 
heat equation for the temporal evolution of distributions 
of temperature. Like the molecular diffusivity and ther- 
mal diffusivity, has its origins in the dynamics and 
interactions that transfer energy. 



B. 



Forces 



k=l 



where: (1) The 5' is a constant of proportionality. Call 
S the conductivity. S is the ability of energy to flow 



The p terms and Uj term in Eq. IjlSfl change by spatial 
movement and are called impelling forces since they cause 
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an energy change. The Nra term, the 



Wj term, 



and 



Tj term in Eq. (|15|) change by temporal differences and 
are cahed nurturing forces since their movement carries 
energy. 



1. Space 

Define a volume V containing Space for the distances 
being considered and throughout which the terms of 
Eq. (|17|l except p terms are constant. The force of Space 
at a point i^sp is defined as 



(18) 



where the negative sign means the force is directed op- 
posite to increasing p. 

Because there is no shear stress in Space, the force ex- 
erted by Space is exerted only normal to surfaces. Con- 
sider a cylinder of cross section ds around a regular, sim- 
ply connected volume with the ends of the cylinder nor- 
mal to Fsp. The ds has a difference of force AFs on each 
end of the cylinder. Allow the height of the cylinder to 
become arbitrarily small. The AFg dFg and 



dF, 



(n • Vp)ds, 



(19) 



where n is the outward unit normal of ds. 

Integrating Eq. (fll^ over a surface, the Space force on 
the surface Fs becomes 



Vp • nds. 



(20) 



The force of Space on a surface is proportional to the 
cross section of the surface perpendicular to Vp. 

The p due to Sources and Sinks can be calculated from 
Eq. H17|l . If all terms of Eq. (|17|l are constant except the 
p, e, and rj terms, the average 77 of each Sink remains con- 
stant over a long period of time, and the average e of each 
Source, each galaxy, remains constant over a long period 
of time, then Eq. 1)17(1 can be solved for a steady-state 
condition. The range of distance from Sources and Sinks 
will be assumed to be sufficient such that the transition 
shock is dissipated and r > rgc- We can consider p at a 
point to be the superposition of the Space effects of each 
Source and Sink. Therefore, for p from the j^^ galaxy 
and fc*^ Sink at the i**^ point (pi), Eq. ((T7jl becomes 



dt 



{r^kP^) = D"^ {r,kP^ ) 



(21a) 
(21b) 



in spherical coordinates, where Vij is the distance from 
the j**^ galaxy Source and is the distance from the k^^ 
Sink to the i'^ point, and pi depends entirely on r.y and 

Tik- 



The boundary conditions are: 

Pi{i'ijT^) — afunctionof distance, only, (22a) 

Pi{rij,t) 0, asr-^cx3, (22b) 

Pi{Tik,0) — afunctionof distance, only, (22c) 

Pi{rik,t) 0, asr^oo. (22d) 

This is analogous to the problem of heat flow in a 
spherical solid with a continuous Source of strength 
at Tij = or with a continuous Sink of strength rjk at 
fkk — 0. By superposition the total pi{x,y,z) may be 
calculated for r^j and r^fc outside the VmP volume 
of the Sources and Sinks as 



P^ir) 



At r,- 



' all galaxies 

\ ^ e 



47rL>2 



all Sinks i i 



(23) 



of the j galaxy , the effect of the other 



Sources and Sinks is very small. Therefore, 



Vp, = 



(24) 



where Vscj is the Tsc for the j galaxy. 

At the j"^ Source, the nearly instantaneous transition 
of the hods and Space into 3D creates a shock wave zone 



of radius 



SWJ , 



where the transition forced Vp is greater 



than VmP- The assumption of constant D is no longer 
valid. In the region Vswj < Tij < rgcj, can be consid- 
ered to be a time average of the nearly digital transitions 
and Vpi = VmP- 

At the Source of the j^^ galaxy — 0, pi{rij = 0) 
is a maximum value pmaxj and Vpi and V^p^ are dis- 
continuous. For periods of the time required for Space 
to flow to Tscj, Cj can be considered a constant. All the 
Space from the Source is flowing through the surface of 
the sphere of radius rscj- Since the volume of a sphere 
cx and e is a volume per unit time, e oc r. Further, this 
proportionality holds for any radius determined by e. If 
ej increases, rgcj increases proportionally as suggested by 
Eq. 121 • Therefore, 



|Vmp| 



(25) 



where Ksc is the proportionality constant which depends 
on how D varies in the VmP zone around the Source and 

Pmaxj ^sc - 

The force Fts exerted on an assembly of hods by the 
Space force due to Sources and Sinks is calculated by 
combining Eq. (|20|l and (|23() and integrating over the 
cross section of the assembly yields 



' all galaxies all Sinks , 



j = l V k=l / 

(26) 

where mg is the Space effective cross section of the as- 
sembly of hods and Gg = (47r)~^. 
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2. Interaction of Space and hods 



where t 



iz and w 



Therefore, 



The hod in 2D is transformed into 3D because of the 
energy pressure of surrounding hods. The transition re- 
leases this energy into 3D. The r acts on Space uniformly 
across the surface of the hod and exerts an attractive 
pressure on Space normal to the surface of the hod. Since 
the hods in 2D are closed loops in a plane, the hods arrive 
in 3D as surfaces oriented flat-to-flat. At the transition, 
Vi = 2dia where di is the initial distance that Space ex- 
tends from a hod's surface. At the radius rh of the hod 
from the Source, Space has expanded and the Pi(ni) < Pi- 
However, within < rh < Tgw, Vp = VmP and ph < Pi 
where ph is the maximum Space density the hod can at- 
tract to its surface. All the Space energy in the shock 
wave zone is within a zone with Vp — VmP- Therefore, 
all the piVi is "bound" to the hods piVi = ra = r^a, 
where Tb is the amount of surface tension used to bind 
Space to the hod. 

Define a Cartesian coordinate system in a volume V 
with the origin at the center of a circular hod with diam- 
eter a, with the z-axis perpendicular to the plane of the 
hod, and with the a;-axis and y-axis in the plane of the 
hod. 

At a radius rh from the Source outside the shock wave 
zone, the maximum distance dh (ho) from the hod surface 
along the z axis within which Vp = VmP is directed away 
from the hod surface has a value 



Ph - Pd 



, Pd = Pi{dh) < Ph, 7-sw < J'h, 



(27) 



where pd is the Space density at dh. 

Hods that are within 2(ih of each other form an assem- 
bly of hods within a volume with a p = pd equipotential 
surface. For Vgw < < Vgc, Vp = VmP everywhere. 

For rh > rsc, volumes more than dh from a hod surface 
have Vp < VmP, P = Pi, Eq- H27|) applies, and T^a < ra. 
Therefore, there is an amount of "free" surface tension 
energy T{ such that 



ra — T\,a -|- Tf a 

Tf = (ih(ph - Pd), T-h > fsc, 

Tb = 2dl - Tf , 

Tf = 0, rh < rsc. 



(28a) 
(28b) 
(28c) 
(28d) 



The Tf is exerted on p^. For the remainder of this Pa- 
per, unless otherwise stated, the region under considera- 
tion is with Pi < Ph. Those hods that are oriented flat- 
to-flat emerge from the Vgc zone as assemblies of hods. 
Other hods become independent. 

Since the energy at the Source upon transition is high 
and since we observe Space to be more extensive than 
the hods, < c?h ^ 1 ho. 

Consider the method of conjugate functions of a com- 
plex variable with a transformation 



cosh u 



X 



sinh u 

^2 



cos^ V sin^ V 



a 

T' 
4 ■ 



(30) 
(31) 



Geometrically, the strip in the w-plane between the lines 
V — and v = tt transforms into the entire t plane. 
The line at it = between t; = and v = n trans- 
forms to the surface of the hod. The lines v ^constant 
are transformed into confocal hyperbolae lines which are 
streamlines of the p.y fleld where Py is the p at the j"^ 
point due to the i*"^ hod. The lines u=constant are trans- 
formed into confocal ellipses with the foci at the ends of 
the hod. Since t acts normally to the surface of the hod, 
u=constant are equipotential lines of pij caused by t [see 
Eq. ©]. 

If the hod is circular around the z-axis, the py stream- 
lines and equipotential lines are also circular around the 
z-axis. The py equipotential lines form oblate spheroids 
with the minor axis along the z-axis. 

Consider a point at a large distance from the hod, r 3> 
a, the u of Eq. (p?n|l is very large. Equation piljl becomes 

+ z^ = d^ where d is the distance from the center of the 
hod. At large d the py equipotential lines are concentric 
spheres. The pij streamlines are the radii of the spheres. 

Since the equipotential volume surrounded by pd is an 
oblate spheroid, the distance d^ (in ho) in the plane of 
the hod to p = Pd is 



dx = 



0.25-0.5, 



(32) 



which is much smaller than a. 

Therefore, the p — equipotential may be approxi- 
mated as a cylinder with height 2(ih and end area = a. 

If Vps > across a hod, the Vp on one side of the hod 
at an inflnitesimal distance Ad from the p = pd surface 
Vpsi(d + dr) differs from the other side Vps2(d — dr). 
Since the Vp is limited to VmP, the force of Space from 
Eq. (|19|l is transferred to the hod surface. Therefore, the 
net force dFg of Space on an element of surface ds of the 
hod is 



dFs 



D^{rfVps2)ds - 
-D^{n • Vps)ds, 



^^(n* Vpsi)ds,(33a) 
(33b) 



t 



cosh(t(;), 



(29) 



where Vps is considered approximately constant over the 
small 2c? distance. This approximation fails in volumes 
close to rsc- 

Unless otherwise stated, when distances greater than 
one ho are implied, d 1 ho and p — ps will be consid- 
ered to be at the hod surface. 

Since the hod is supple, the application of dFs on ds 
of the hod in the x-y plane will cause ds to move nor- 
mal to the hod surface. If dFg — everywhere else, a 
deformation, du, of the hod surface normal to hod sur- 
face will occur. Since a and t are constant, the edge of 
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the hod must change shape. A small deformation allows 
movement in only the z direction. The edge stays in the 
x-y plane. However, since edges of a free, individual hod 
are not constrained and r is very strong, the hod mem- 
brane does not oscillate normal to the hod surface as it 
would for a stretched membrane. Therefore, dFg will 
cause the hod to tend to conform to the Vps stream- 
lines. Since Space has a relaxation time, the movement 
of hods will cause the streamlines to be continually shift- 
ing. A slight delay, caused by l and by 5, in adjusting to 
a change in the streamline will allow interesting effects. 
If Vps • n 7^ for a brief instant on only one edge of 
the hod and if the other edge experiences a force par- 
allel to the surface Vps • n — 0. The hod will travel 
along the streamline at the maximum allowed velocity. 
If the streamline the hod is following has a very sharp 
bend, such as at the edge of another hod where the hod 
is much stronger than anywhere except near the Source, 
the inertia energy of the hod will let it continue. For 
an instant, Vph • ''^ 7^ 0. The hod will slow or even 
reverse direction. This allows the transfer of energy be- 
tween hods. 

If the Vp field of Space varies little over the surface of 
the hod, the deformations become a general shaping of 
the hod causing the hod to tend to conform to the force 
lines, or streamlines, of p. The hod will move to make 
Vp • n — on the surface of the hod. 

If Vp • n = from some direction, the hod will expe- 
rience zero Space impelling force from that direction. If 
the Vp field changes little over a range of several hod di- 
ameters, small deformations with Vp • n ^ will cause 
the hod to have a forward inertia. 



Gravity 



The boundary conditions are: 

Ui{dij,0) — a function of distance, only, (36a) 
Ui{dij,t) Oasr— > 00. (36b) 



Analogus to Section III B II 



(37) 



The \Uij\ is highest in a volume with a given Tf > 



value when the k and j hods have a shared p = 
Pd oblate spheroid surface and are orientated flat-to-flat. 
Since dij is the distance to the center of the hod, rfy ^ 
and a singularity problem is non-existent. 

Chose two assemblies of hods where Nt is the num- 
ber of hods in one assembly and Ns is the number of 
hods in the other assembly. Let each assembly have a 
very small diameter compared to the distance rts =< 
^2^=1 J2k=i ^kj > where < > means "average of". The 
force of gravity of the hods in an assembly is greater on 
each other than the force of gravity of hods in the other 
assembly. Therefore, the hods of an assembly will act as 
a solid body relative to the other assembly of hods. Be- 
cause hods act on Space, Space acts on Space and Space 
acts on hods, only the local T{ for each hod determines 
the N Tf a on the hod. Therefore, the force of gravity 
on each assembly due to the other assembly is 



(NsTk) (iVtTft) rts 



(38) 



where Tfs is the Tf for the assembly with Ns hods and rft 
is the Tf for the assembly with A^t hods. 

For the simplified case of two symmetric distribution of 
hods with a distance ri2 between them and with Tfs = 



Call Fgk the force of gravity on the k hod 



S ^' 



(34) 



where the negative sign means the force is directed oppo- 
site to increasing U and the sum is over the other (j'^) 
hods. 

Define a volume V containing Space for the distances 
being considered and throughout which the terms of 
Eq. H17() except the U terms are constant. By the Princi- 
ple of Superposition, this allows the examination of just 
one force. Equation ifTTjl can be solved for a steady-state 
condition. We can consider p at a point to be the super- 
position of the effects of each hod. Therefore, Eq. H17|l 
at the i"^ point due to the j"^ hod becomes 



^^id,,U,) = §V^(d,,C/, 
in spherical coordinates. 



(35) 



mgimg2 
-Lf o ri2. 



(39) 



where G = (47r)~^ (Kf^a)^, mgi — NtTn, and mg2 = 
As Tfs. The experimental determination of mgi and mg2 
must use a method depending on Uij and no other factor 
such as motion which will involve La forces. 

Note, in volumes within a sphere with a radius of rsc 
from a galaxy's Source, T[ — and, therefore, Fg — 
and TOg = 0. In volumes with rsc < ^'h, Tf increases as 
increases. 



C. Particles 

Define a particle to be an assembly of hods surrounded 
by a p = pd equipotential surface (ES) that forms a con- 
nected volume around the hods in the particle within 
which Space is totally bound to the hods by Tb. There- 
fore, Ps outside the ES equals pd within the ES of the 
particle. If the ES is simply connected, the particle is 
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called a simple particle and the distance for gravity po- 
tential calculations is to the center of the particle. The 
potential on a surface element of the ES is determined 
by the distribution of hods inside the surface. However, 
the volume occupied by the bound Space is the same for 
each hod. Therefore, the hods are uniformly distributed 
over the volume inside a simply connected ES. Since the 
hods of a simple particle have a minimum distance be- 
tween them, the hods are tightly bound and respond as 
a single entity. By the Principle of Minimum Potential 
Energy, the hods of a simple particle have a more stable 
structure with a maximum Nr^^a per unit surface area 
ratio (P/S). 

Hods are unable to join outside the Source region to 
form a particle but may be able to join an existing par- 
ticle. The existing particle must have sufficient 7V(rfQ;) 
to overcome the ia of the joining hod. Also, this implies 
there is a minimum number A^min of hods which must 
be combined before other hods may join it outside the 
Source region. 

When the distance between the flat of the hods < 2rfh, 
the Vp between them will change sign and become VmP 
which will cause a net repulsive force between hods. 
Therefore, the distance between the flat of the hods will 
stabilize at a distance 2dh when the repulsion due to VmP 
will equal F^. If the distance between the flat of the hods 
> 2dh, the hods will be independent. 

Also, since rfh 1 ho, hods may form in a column 
arrangement. Call this type of structure Family 1 (Fl) 
Type 1 (Tl) structures with the length L — 2Acdh where 
A"c is the number of hods in the FlTl structure. The po- 
tential can become large enough that hods can be added 
without being in a Source r^c zone. The P/S of the FlTl 
structure will be least when 2A'niinC?h = 1 ho The addi- 
tion of hods to the column also produces stable, cylinder 
structures since each additional hod is oriented flat-to- 
flat. Therefore, the addition of hods to column length 
decreases the P/S. However, for a large number of hods, 
the outer hods are easier to remove. The direction paral- 
lel to the surface has Vp»n = so Fg = and maximal 
speeds can be achieved in a direction parallel to the sur- 
face. This structure is a photon. The direction along the 
axis has a one-hod surface. Hence, Fs along the cylinder 
axis is as if only one hod were present. 

If the inertia and external forces cause the hods in the 
FlTl to start to loose overlap, the sharp Vp fleld at the 
edges of hods will push them back. The total overlap 
condition is a minimum potential energy position. 

Any deformation in the surface from a ripple in Space 
will tend to bend the hod surfaces to lose total overlap. 
Since this causes a force to reestablish overlap, this struc- 
ture is much more rigid than a single hod. Because this 
structure is rigid, a Vp across the structure can impart 
a torque on the particle. 

As more hods are added, the Nc{Ta) increases. The 
energy of the photon hv is 

hi/ = Afc(™). (40) 



Ripples in Space will cause ripples in the surface of the 
outer hod. This will slow the photon structure and cause 
a shearing F^ on the outer hod. Thus, if A'c > Amin in 
a photon, the photon can loose hods or gain hods. 

Two cylindrical photons may join as side- by-side cylin- 
ders caused by Fg. Parallel to the surface, Vps • n — Q 
for both. 

Also, two cylindrical photons may join in a T formation 
if L > 1 ho as the Principle of Change requires. By 
extension, if L « 2 ho, there may be other cylinders 
added to form II, HI, etc. structures. This requires the 
photons in the vertical columns to have approximately, 
if not exactly, the same number of hods. Otherwise, one 
or more of the photons will experience a Vp • n ^ Q. 
There will still be a direction where Vp • n = 0. Call 
this structure a Fl, Type 2 (T2) structure. Perhaps this 
is the structure of a lepton. 

Consider a F1T2 structure made of several columns 
and a cross member at each end. If an end column does 
not overlap the cross member column, an extremely un- 
stable condition exists. The end column would be freed. 
If a cross member is longer than the width of the row 
of columns, the particle could acquire another column. 
The most stable structure would be when L is an integer 
multiplier / of the average diameter of a ho 

L = 2/(-)i/2. (41) 

TT 

Similarly, a third photon can join the T structure so 
that the three are mutually orthogonal, as the Principle 
of Repetition allows. However, Vp • n ^ in all direc- 
tions. Call this structure a Fl, Type 3 (T3) structure. 
A cross cylinder on a cross cylinder of a F1T2 structure 
makes the FITS box structure very dense and very stable 
since there is little "free" surface (unattached hod faces) 
compared to the rta energy contained and P/S is mini- 
mal. Perhaps this is the structure of quarks, quark stars, 
and black holes. 

In small structures, the F1T3 cube has the least P/S. 
However, this allows "corners". In larger structures, a 
rounding effect such as found in crystal structures pro- 
duces a more spherical ES which will have a lower P/S. 
This could be a mechanism to fracture some combina- 
tions (particles) and cause faceting in all combinations 
like in crystals. 

Another possible structure of hods with a position 
of least P/S and stable Vp fleld is with each inserted 
through and normal to each other (see Fig. [SJ. Call the 
particles formed with this structure Family 2 (F2) par- 
ticles. The particle so formed has a direction along the 
line where they join where Vp»n = can occur. These 
may be the Muon family of particles. 

By a similar argument, with less probability, three 
hods may come to form a particle with three mutually 
orthogonal hods (see Fig. ^ . Call the particles formed 
with this structure Family 3 (F3) particles. These may 
be the Tau family of particles. 

However, a 4**^ mutually orthogonal hod construction 
in 3D outside the immediate area of a Source is not fea- 
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FIG. 2: Diagram showing the structure of the Family 2 
(Muon family) fundamental particles. 




FIG. 3: Diagram showing the structure of the Family 3 (Tau 
family) fundamental particles. 

sible since one of the hod's surfaces will face another hod 
surface. 

Therefore, because the hod defines 3D, there are three, 
and only three, famihes of particles. 

In our volume of the universe, particles of the same 
type have consistent mass. Particles become larger from 
gravitational attraction of additional hods as the Princi- 
ple of Repetition allows. The observed consistent mass of 
like particles implies the existence of the Principle of Neg- 
ative Feedback to limit the size. For T2 and T3 particles, 
if L > /(a/7r)^/^ from Eq. attracted hods may be 

added to the particle. As noted previously, the attracted 
hod will be traveling at the maximum speed which im- 
plies a maximum T . Thus, the attracted hod has energy 



great enough to break the gravitational bond and facet 
the outlying hods from the particle. If L < /{a/7r)^/^, 
the energy of the attracted hod may cause the column to 
facet. 

In a volume with slow changing Pi, Pd = Pi- Therefore, 
L is a function of pi. In a volume with rapidly chang- 
ing Pi, Pd may lag the change in pi because rb may be 
absorbing or releasing Space. Such may be the situation 
for neutron stars, possible quark stars, and black holes, 
which exist near the Source of the Galaxy. 

If the particles are to have more hods (energy/mass) 
as the Anthropic Principle requires, there must be an- 
other binding mechanism less sensitive to faceting as the 
Principle of Change requires. The flat side of the end 
hods on a photon would be attracted to and could be 
bound to other particles. If quarks are the largest Fam- 
ily "type" particles, baryons must be quarks and leptons 
joined by another means such as directly to each other or 
by photons. The particle will be in a connected volume, 
but not a simply connected volume. By the Principle of 
Repetition, the atom is electrons, neutrons, and protons 
joined by photons. The photon is similar to a rod in these 
structures. 



D. Sink characteristics 

If a Source ceases to erupt or if some matter is repelled 
to the low Vp distances between Sources, gravity will 
cause the hods to coalesce. In a volume where Vp be- 
comes near zero, the attraction of hods will cause a high 
hod density. If the number of hods is greater than in a 
black hole in this intergalactic volume, the intense force 
on the supple hods may cause the hods to form spheres 
like bubbles. With increased pressure from a huge num- 
ber of hods, the hod can go only out of our universe and 
turn into a four dimensional 4D object like 2D turned to 
3D. 

Since the Sink requires mass to form, the Sink's age 
is considerably less than galaxies. Therefore, the delay 
before matter and Space is sent into the 4D allows 3D to 
accumulate mass. 

The Sink requires a minimum mass to send hod bub- 
bles into 4D. Therefore, the mass Mgink around the Sink 
has an upper limit. If the Sink is sending hods to 4D, 
then a negative feedback condition is formed and the rate 
rj of transition of the hods is dependant on the rate of ac- 
cretion of hods. Therefore, ry oc Mgin'k- 

The Vp acts in the same direction as the gravity force. 
The Sink can send hods to 4D if the hods are close. 
Therefore, the Vp = Vmp around the Sink like around 
the Source. 
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FIG. 4: Diagram depicting the regions of a galaxy in relation 
to a hypothetical rotation curve. The lines in the Saturation 
Region depict the different shapes possible. 



III. GALAXY MODEL 
A. General 

For the purpose herein, a galaxy's RC is divided into 
five regions [23 as the r from the Source increases (di- 
agramed in Fig. a center region (CR) [s^, l4H| . 
a Keplerian region (KR) Il6| . a rising region (RR) 
which is in the outer bulge region, a saturation region 
(SR) which is in the disk region, and a detached region 
(DR) which is outside the galaxy (herein, collectively re- 
ferred to as the "lettered regions"). The distinct nature 
of each region of the RC implies different physics exists 
among the regions. Therefore, some of the parameters 
of the equation of motion have suddenly become more or 
less significant. Each of the lettered regions is separated 
by a transition region in which the physics is changing 
from one region to the next. The CR and KR are sep- 
arated by a transition region Tc-k E 111 IM HI • The 
KR and RR are separated by a transition region Tk-r 
in which the RC has a slight decline followed by a rapid 
rise. A similar transition region Ti._s appears between 
the RR and SR. Particles in DR are not in orbit around 
the galaxy center. 

Dividing into regions with different, predominant 
physics has been done for models of other phenomena. 
For instance, the physics of semiconductors uses such a 
model to describe differing current-voltage characteristics 
of the differing regions of transistor operation [sif . 

The T[ — out to Tsw and only unbound hods and 
photons exist near the Source. By Eq. (|25|l the size of 
e is proportional to Vgc and, therefore, to the diameter 
of the Source region. As the hods and photons travel 
outward beyond rgc, Tf > and particles form. 



As the particles' radius r from the Source increases, 
larger particles can form. Outward from the Source, pho- 
tons coalesce to form larger particles. In the Big Bang 
model, the physics of the formation of particles is out- 
ward in time. In the CUM, the physics of formation 
of subatomic particles is outward in distance from the 
Source. Eventually, hydrogen forms and is pushed out- 
ward by the Space force. As hydrogen coalesces, larger 
atoms are formed which fall back to the center as the 
gravitational force becomes larger than the Space force. 
Gravity attracts the densest particles into the center. As 
the particle approaches rgc, the declining Tf and L causes 
the T2 and T3 particles to disintegrate to photons and 
re-radiate outward. 

Consider a AV containing a test particle with N hods 
such as a star or hydrogen gas in a galaxy. The iner- 
tial mass is Nia, the gravitational mass rUg is Nrta, 
and the cross section mass nis is the effective cross sec- 
tion subject to the Space force of the hods in AV. For 
simplicity, the number of hods in the particle is assumed 
to be constant over time. Since stars are changing the 
number of hods by radiation emission and are changing 
the ms/rrii ratio through changing elemental composi- 
tion, this assumption is most nearly true when the AV 
contains hydrogen gas. Therefore, the Hi RC is preferred, 
where possible, to reflect the forces infiuencing v. 

Assume the particle in AV originated in the galaxy 
(db = 0) and the particle in AV is in equilibrium [E^ = 
NTfo). The coordinate system center is placed at the 
Source and is aligned in a constant orientation relative to 
our view. The velocity and rotation of the galaxy relative 
to the universe is not part of v. Combining Eqs. I|23|l . 
(|S7|l . ©, and ifT^ into Eq. ((T3|l . rearranging terms, and 
considering only radial terms yields 



rr - 



mg GM nis Gge 



TO,, 



TO, 



r • ao. 



(42) 



where the double dots above the r mean the second 
derivative with respect to time, r is the radial acceler- 
ation of the test particle, G is the gravitational constant, 
Gs is the cross section constant, M is the mass of the 
galaxy inside a sphere with a center at the Source and a 
radius of r, e is the Source strength of the galaxy. 



other galaxies 

E 



TTlgGMj 



TTlsGs^j 



E^'""^^ rrigGMk + msGs\rik\ 
■ 



N 

i=l 



(43a) 
(43b) 



th 



Mj is the total mass of the j galaxy, is the Source 
strength of the j"^ galaxy, is the distance from the i 
test particle to the j"^ galaxy's Source, Mk is the total 
mass of the fc"^ Sink, 77^ is the Sink strength of the k^^ 
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Sink, m, is the inertial mass of the test particle, and Tik 
is the distance from the z*'' test particle to the k^^ Sink. 

If Co ^ in a galaxy, then the center of mass may be 
slightly displaced from the Source. The v is measured 
along the major axis from our view. Therefore, the r is 
directed along the major axis. The form of Eq. H43|l for 
baryons is of a Newtonian force with an effective mass 
less than the total mass in a galaxy. Also, the force on 
FlTl (photons) particles passing by the galaxy differs 
because of the orientation of the particle relative to the 
galaxy. 

In our part of the Galaxy, = mg to within one part 
in 10^^ |52l |. However, the measurement of rrig depends 
on i^g and, therefore, rf that depends on p. The ratio of 
the diameter of the earth divided by the distance from 
the center of the Galaxy suggests much greater precision 
is required for earth bound measurements. Therefore, 
the TOg/m,, ratio will be explicitly stated in this Paper 
where applicable. 

The notation convention used herein is as follows. The 
italic, normal size letter will denote the parameter, as 
usual. The upper case M means the total mass inside 
the radius denoted by the sub letters. The sub letters 
denote the radius range by region identification. The 
use of "max" or "min" in the sub letters will denote the 
maximum value or minimum value the parameter has 
in the region, respectively. For example, the radius 
(kpc) in the KR which is between the radius rtckmax at 
the maximum extent of the Tc-k region and the radius 
^kmax at the maximum extent of the KR region. Hence, 
T-tckmax < ''k < J'kmax- Thc notation Mkmax mcaus the 
total mass inside the sphere with radius rkmax- This is 
the total mass of the CR, Tc_k, and KR regions. 

The CUM suggests the change in physics among the 
lettered regions derive from a change in the particle 
TOs/mi, in the forces (i^ts and -Ftg), and in r»ao. The let- 
tered regions change is shared by all galaxies and, hence, 
may be used to derive relations among galaxies. 

The following sections discuss each region. Some sec- 
tions have data that can be compared to the equations 
developed. 



B. CR to KR 

The CUM posits that the radius of the CR is less than 
a few light hours around a Source at the center of a galaxy 
and is less than rsc- At a radius of less than rg^, there 
is no V and the problem of a singularity in Eq. H42|l at 
r = is nonexistent. Since Tf decreases as the r decreases 
below Tsc, the gravitational term is decreasing and the e 
term dominates to repel particles. 

The unbound hods and photons in the CR are repelled 
away from the Source. The e/r close to the Source is 
large. Therefore, only exceptionally dense particles with 
very small ms/m^ are closest to the Source. Black holes 
which may be F1T3 particles and other very dense matter 
fall into the Tc_k of galaxies. All other matter is repelled 



outward. Black holes form a shell around the CR 29ll41j. 

Black holes may increase nis/fn^ by collision and 
merger or may loose angular momentum such that they 
fall into the r < r^c zone. The result of the L ^ and 
Tf — > is a multitude of X-ray point sources without ac- 
companying large variation at longer wavelengths near 
the center of the Milky Way. Note the force of gravity 
must be reduced to allow the photons to separate in a 
burst. The diffuse component of the X-ray sources may 
be other low ms/wt particles such as iron nuclei or neu- 
tron stars with F1T3 centers and references therein]. 
At and near r^c the brightness of the radiation will in- 
crease relative to star luminosity due to the reclaiming 
of the dense particles. 

If mass does fall into the Tc-k, the higher mass 
A/tckmax in the Tc_k, causes more infall into the CR. 
The resulting increase in re-radiated photons causes a de- 
crease in Mtckmax- From Eq. (|42|l . the e determines the 
stable value of Mtckmax- This is the observed feedback 
controlled mechanism required by Merritt and Ferrarese 
[35I I as proposed by the CUM. 



KR 



The start of the KR is caused by a sudden and rapid 
change in particle type. Almost all the very dense parti- 
cles ( low TOs/r/it) are in the Tc-k shell. As rk increases, 
the number of dense particles per unit radius decreases, 
the resident particles's rris/m^ term is greater, and the 
Aftckmax is much larger than the gravitational effects of 
particles in the KR upon each other. The observed Kep- 
lerian motion of particles in the KR 15, 16] implies the 
mass Mk within a sphere of radius rk is approximately 
-Mtckmax- Also, siucc Fg oc Tf of the test particle, the 
observed Keplerian motion of particles in the KR implies 
the Tf of the test particle is nearly constant in thc orbit 
of each test particle. 

The TOs/mt term is causing lighter and denser parti- 
cles to move out of the KR. The ms/m^ for hydrogen 
and the lighter elements is much higher than the massive 
particles and the mg is lower. Thus, for hydrogen in the 
KR the predominant term is the e term which pushes 
hydrogen and lighter particles out of the KR. The only 
lighter elements in the KR are gravitationally bound to 
the denser particles. Black holes are either in or on their 
way to Tc_k- Therefore, Eq. H42|l for particles in the KR 
becomes 



2 .. . TOg GMcff 

I'k = '"kT'k H Ak, 

Tk 



where 



Mtcki, 



G 



(44) 



(45) 



Wk is the V for particles in the KR, and K]^ is a relatively 
small constant accounting for the term. An upper case 
K followed by a sub number or letter denotes a constant. 
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The Tk in Eq- (|44|l has a value to cause elhptical orbits 
for the particles which stay in the KR. 

As rk increases, ^ in a Keplerian decline. Since 
differing particles have different ms/m^ ratios such as 
stellar material and hydrogen, the minimum Wk for dif- 
ferent particle species is reached at different rk. 

The sign of rk • flo is opposite on opposite sides of the 
galaxy. Also, the distance between rk on one side to rk 
on the other side is of the order of a few kpc. Compared 
to the other terms, rk • is small. Therefore, the resid- 
ual of the rk • Qo from each side of the galaxy is very 
small. Therefore, the luminosity curve as a function of 
rk will be only slightly different on each side of the cen- 
ter. Therefore, averaging the rotation velocity from side 
to side reduces the effect of relatively small rk • flo term. 

Since the Mtckmax is concentrated in a small volume, 
the KR is spherical. Galaxies are usually inclined relative 
to our view. Because KR is spherical, the diameter of the 
KR is measured as the distance along the major axis of 
the ellipse (from our view). 

Since the ms/m^ particle type differs with radius, this 
Paper posits that in the high force volume of the KR 
the ms/m^ particle types are in strata and the strata's 
radius depends on e. If each mg/wt particle type has a 
different mass to luminosity ratio, the surface brightness 
/ versus r curve {I — r curve) along the major axis has 
discontinuities or inflection points. Define r^d (in pc) as 
the average radius along the I — r curve in each direction 
from the Source to the first measured discontinuity or 
inflection point. Therefore, Eq. H25I) suggests r^c oc r^d 
and 



(46) 



where iiTed is the proportionality constant. 

Further, the abrupt change in ms/m^ particle type oc- 
curs in all galaxies. Therefore, r^d may be used to com- 
pare parameters among galaxies. 

Evaluating Eqs. and at r^d, solving for r^d 
with the net effect rk • flo ~ when values from each 
side of the galaxy are averaged yields 



mgGM{r,d) 



(47) 



where M{r^d) is the mass within a sphere of radius r^d 
and w(red) is the rotation velocity of particles at r^d. 

The term in the braces is the difference between the 
measured t;^(red) and the Newtonian expectation. 

Substituting Eq. (gSl into Eq. (gSJ yields 



G m 

Moff = Mtckmax - -T^Ked—red + Ki^B, (48) 



where Ki^s is the intercept of the M^s — r^d plot. Theo- 
retically, i^ioff = and is a result of measurement error. 



1. M, — Ted data 

The kinematic measure of Af, uses the relation ra^/G. 
Therefore, from Eq. (|44ll M, « Mgg, where the mass in 
the KR is relatively much smaller than the mass in the 
Tc-k- 

To test Eq. (gHJ, M, data using kinematic methods 
and distance data for 20 sample galaxies were taken from 
Merritt and Ferrarese js^. The data for the 20 sam- 
ple galaxies are listed in Table The data points are 
the measured values for the galaxies for which the SBH 
sphere of influence (rh in Merritt and Ferrarese [s^) has 
been resolved. The morphology data were obtained from 
the NED database The sample includes E, SO, and 
S type galaxies. 

The red data were calculated from the distance in 
Refs. [35I I and from the average angular (observed) radius 
ra (arcsec) . The r^, for each sample galaxy was measured 
from Hubble Space Telescope (HST) images listed in 
Table m A Visual Basic program (see Appendix ^ was 
written to extract the I — r curve and calculate ra. The 
selection of each image was based upon being in the V 
band, upon having a sufficient area of the center of the 
galaxy in the image so ra on both sides of the galaxy 
could be measured, and with sufhcient exposure time to 
have the / value [80| at ra sufficiently large. The I — r 
curve obtained from the HST images using the WFPC2 
instrument are shown in Fig. El NGC 4607 which is listed 
in Merritt and Ferrarese [33 lacked a satisfactory image. 
Since the goal of the program is to select a pixel at a 
discontinuity, image analysis techniques such as neigh- 
bor pixel averaging, azimuthal averaging, and combining 
more than one image into a single data set were avoided. 
The Visual Basic program has alternate means to deal 
with peak (center) saturation, bad high pixel values, and 
bad low pixel values. 

The straight line in each I — r curve in Fig.|Slis drawn 
between the pixels the Visual Basic program selected to 
calculate ra. Note, the line is at differing / levels among 
the galaxies and the slope of the line indicates asymmetry 
in / value for the calculation of ra in the same image. 
Therefore, the selection of ra is not based on selecting an 
isophotal value. The ra is not an isophotal based effective 
radius. 

Figures El shows plots of M, versus red of the selected 
galaxies. The straight lines were drawn by (1) identi- 
fying the half of the galaxies furthest from the origin, 
(2) drawing lines from the origin to the identified galaxy 
points, (3) noting six lines could be matched to the iden- 
tified galaxies, (4) the galaxies were organized into clas- 
sifications according to their proximity on the graph to 
a line, (5) finding the slope and intercept using the least 
squares method for each classification group of galaxies. 
The differing slopes of the lines marked by filled dia- 
monds, filled squares, filled triangle, "X", star, and 
suggest either differing ms/irig values, differing K^d val- 
ues, or differing Mtckmax in Eq. among the galaxy 
classifications. That the data points are on a straight line 
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FIG. 5: Plots of the J — r profile for the sample galaxies. The straight line is drawn between the chosen pixels. The slope on 
the line indicates the asymmetry. 
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TABLE I: Data of the galaxies in Fig.^l 



Galaxy" 


Morphology 


HST 




Distance' 






line^ 






image 


filter 


Mpc 


pc 


10' 




n 


11459 


E3; LINER 


U2BM0102T 


F555W 


30.3 


99.5±14.7 


4.6 


± 2.8 


2 


N0221 


cE2 


U2E20309T 


F555W 


0.8 


3.7± 0.4 


0.039± 0.009 


1 


N2787 


SB( r)0+ LINER 


U39S0304R 


F555W 


7.5 


14.6± 3.6 


0.41 


+0.04 
-0.05 


2 


N3031 


SA(s)ab: LINER Syl 


U32L0105T 


F547M 


3.9 


5.5± 1.9 


0.68 


+0.07 
-0.13 


3 


N3115 


SO- 


U2J20B04T 


F555W 


9.8 


49. 1± 4.8 


9.2 


± 3.0 


4 


N3245 


SA(r)-0- 


U3MJ1103R 


F547M 


20.9 


35.5±10.1 


2.1 


± 0.5 


3 


N3379 


El: LINER 


U2J20F03T 


F555W 


10.8 


51.7± 5.2 


1.35 


± 0.73 


2 


N3608 


E2 LINER 


U2BM0502T 


F555W 


23.6 


56.5±11.4 


1.1 


+ 1.4 

-0.3 


2 


N4258 


SAB(s)bc;LINER Syl.9 


U6712101R 


F547M 


7.2 


12. 1± 3.5 


0.39 


± 0.034 


2 


N4261 


E2-3: LINER 


U2I50204T 


F547M 


33.0 


189.4±16.0 


5.4 


± 1.2 


2 


N4342 


SO- 


U32Q0206T 


F555W 


16.7 


36.9± 8.1 


3.3 


+ 1.9 

-1.1 


3 


N4374 


El: LERG LINER 


U34K0102T 


F547M 


18.7 


240. 2± 9.1 


17 


+ 12 
-6.7 


3 


N4473 


E5 


U3071503T 


F555W 


16.1 


49. 1± 7.8 


0.8 


+ 1.0 

-0.4 


1 


N4486 


cEO 


U2900104T 


F547M 


16.7 


32. 1± 8.1 


35.7 


±10.2 


6 


N4564 


E6 


U3CM7202R 


F702W 


14.9 


26.0± 7.2 


0.57 


+0.13 
-0.17 


2 


N4649 


E2 


U2QO0301T 


F555W 


17.3 


52. 7± 8.4 


20.6 


+5.2 
-10.2 


5 


N5128 


SOpec 


U4100106M 


F555W 


4.2 


33. 2± 2.0 


2.4 


+3.6 
-1.7 


3 


N5845 


E 


U3070903T 


F555W 


28.5 


167.3±13.8 


2.9 


+ 1.7 
-2.7 


2 


N6251 


E:LERG Sy2 


U2PQ0701T 


F555W 


104 


206.9±50.4 


5.9 


± 2.0 


2 


N7052 


E 


U2P90108T 


F547M 


66.1 


400.3±32.1 


3.7 


+2.6 
-1.5 


1 



""I" means IC and "N" means NGC. 

''The distance is as stated in the reference for M, . 

•^The r^d (kpc). The error is ± 0.1 arcsec. 

"^The M, was obtained from Merritt and Ferrarese l35l . 
=The line is the "n" listed in Table |3 
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FIG. 6: Plot of the measured mass of the theorized supermas- 
sive black hole M, (10* Mq) at the center of galaxies versus 
the distance along the major axis of a HST image to the first 
discontinuity or inflection point r^d (pc) using the distance 
stated in Merritt and Ferrarese Issl . The filled diamonds, 
filled squares, filled triangle, filled circles, open diamond, and 
open square are data points along a linear M, - r^d line using 
data from Table Q and correspond to n = 1, n — 2, n — S, 
n = 4, n — 5, and n = 6 as listed in Table HTl respectively. 
The error bars (see Table HJ were omitted for clarity. 



implies the m-s/mg ratio (particle type), i^ed, or Mtckmax 
is discrete for an e value. Since the lines of n=l to n=3 
(the lines with more than one point) nearly intersect at 
(r,d,M.)=(0,0) and Mtckmax > 0, Eq. gHJ implies 

Mtckmax = K,^e = K.^K^ATed, (49) 

where K^e is a proportionality constant. The K]^ of 
Eq. H44|l produces a small offset in the intercept of the 
lines of Eq. 
Therefore, 

M. « Mcff = -FCsoftr.d + i^icff, (50) 

where 

K,,s = K„K,d - ^K,d— (51) 

G TOg 

and KscS is the slope of the M, — r^d plot. 

For a spherically symmetric particle, nis increases with 
the square of the radius and m, increases with the cube 
of the radius. For a given density, larger particles have 
lower TTis/m,,. Hydrogen is found in the outer regions of 
a galaxy and heavier elements arc found in the inner re- 
gions. Stars larger than our sun are found both inward 
and outward of us. Therefore, hydrogen has a higher 
ms/rrii than iron. Black holes have a much lower ms/rrii 
than other stellar particles. Therefore, the largest cross- 
section area on which Fg acts is on the area of nuclei 





TABLE 11: 


Data of the lines in Figures 


m 


n" 


hue 


Corr." 




7 

A'icff " 








10** Mq kpC-1 


lU M0 


1 


filled diamonds 


0.99-f 


9.0± 0.3 


nil 1 A r\'7 

0.11±0.07 


2 


filled squares 


0.92 


27 ± 6 


0.07±0.8 


3 


filled triangles 


0.99+ 


70 ± 2 


0.1 ±0.3 


4 


filled circle 




187 





5 


open diamond 




391 





6 


open square 




1113 






"The integer denoting the place of the line in the order of increas- 
ing ifscff • 

''The identification of the line symbol in Figures [51 
■^The correlation coefficient. 

''The least squares fit slope of the M. — r^d lines. 
■^The least squares fit intercept of the M. — r^d lines. 
■^AU lines are calculated as the least squares best fit for points and 
(r,d,M.) = (0,0). 

rather than on atoms, compounds, or larger assemblies 
of matter. Gravitationally bound assemblies of matter of 
differing sizes and the same relative elemental composi- 
tion will experience the same Meff. Also, assemblies of 
matter with differing relative elemental composition will 
have differing RCs. Differing rris/mi^ cause differing at 
a given r. For instance, the Hq, (found near hot stars) 
and Hi RCs are different at small r because the particles 
are very different. However, in the disk as r increases, the 
stellar elemental composition is of successively lighter el- 
ements. Therefore, the H^ and Hi RCs become similar 
at larger r. 

The slopes among the lines in Table|n]obey the relation 

logio(^) =A;„n + ifi„, (52) 

where S'lope is the slope of the lines from the parameter 
relationship; Unit is the units of measure of S'lope; and Ksn 
and Kin are the slope and intercept of the linear relation 
of Eq. (|52(l . respectively. For the M, — r^d relationship, 
^lopc = K,,s, Unit is lO^Mo kpc~i, K,n - 0.41 ± 0.01, 
and Kin = 0.58 ± 0.04 at one standard deviation (1 a). 
The n is an integer of one to six depending on the line 
from Table ITU The correlation coefficient and F test [sj 
for Eq. lO are 0.999 and 0.992, respectively. Note, the 
lowest slope corresponds to n=l of Eq. H52() and the lines 
of Figures ini intersect at a point. Note, S'lope is nearly in- 
dependent of the distance used to calculate the M, and 
Ted if the same distance is used for both. The values of 
logio(e)« 0.434 is within the 2.3 a of log]^Q(i4rsn)- There- 
fore, restating Eq. (|52|l as a strong Principle of Repetition 
yields 

X,cff = ifsce", (53) 

where K^c = (3.7tEJ;J|) x lO^Mekpc'^ at 1 a. 

Combining Eqs. H50|l and H53|) and setting K-^ff = 
yields 

M. =if,ce"r,d. (54) 
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The standard deviation of the function are the er- 
ror bars of M, Usted in Table |T| The = 17 and the 
probabiUty value of the Px ~ 60%. The data does 
not invalidate Eq. (|54|l at the significant 5% level. The 
and P-^ were calculated using only the galaxy data 
(not the {r^d,M,) = (0,0) point). The Af, is a function 
of two parameters, r^d and n. The r^^ is a parameter 
of each galaxy, the n is a galaxy classification parameter, 
and the e implies the Principle of Repetition is applicable 
to the M, parameter. Given a set of parameters (n, r^^), 
a unique value of M, can be calculated. 

A feature of Fig. is the lack of data points in the 
upper right of the plot (galaxies with relatively high red 
and high M, values). 



D. RR 

In the KR as increases to its maximum value 
(''kmax), the rotation velocity Wkmin decreases to its min- 
imum value and the Space and gravitational forces bal- 
ance. This implies 

Mefi ^ 0, (55a) 

^kmin ^ ''kmax(^kmax ~ ^kmax • 0,o) ■ (55b) 

Therefore, in Tk_r the motion of the test particle is pri- 
marily governed by the force from external galaxies. Par- 
ticles such as hydrogen are being ejected outward from 
the KR with a large radial velocity can cross the Tk-r- 
Particles such as dense stars in nearly circular orbits on a 
slow inward journey may lack the necessary radial veloc- 
ity to cross the Tk_r. The Tk_r is a centrifugal barrier to 
inward movement of particles. Therefore, external galax- 
ies can have a large influence on the mass distribution and 
total mass in a galaxy. 

If Wkmin = in a galaxy, then the CR and KR are de- 
coupled from the RR and SR. The w at r > rkmax is as if 
there were no mass inside rkmax and, for non-relativistic 
calculations, the particles in the RR and SR of an in- 
trinsic galaxy [a^ — 0) are subject to only Newtonian 
forces. 

Two changes in the physics occur to cause the RR. The 
first change in the first RR (FRR) is to change the par- 
ticle elemental species which changes mg/mt while main- 
taining a spherical region shape. The second change in 
the second RR (SRR) is to change the dr increment shape 
from a spherical shell to a cylindrical shell while continu- 
ing the changing mg/m, . The resulting equations depend 
on the V reaching a minimum value ( dw /dr = 0) of the 
RC in the KR for the particle species examined. 

That shapes of astronomical particle assemblies are ei- 
ther spheroidal or disk is well established. This Paper 
takes the disk shape as having cylindrical symmetry. The 
well-known peanut shaped profile of a galaxy bulge is a 
varing height, cylindrical symmetry with a spheroid core. 



1. First RR 

The rotation velocity Wr of a particle in the RR is 
greater than Wkmin- If the radius rtkr of a particle's orbit 
in the Tk-r decreases to less than rkmax, the attractive 
force decreases because M is decreasing. Thus, the ef- 
fect of e which is spherically symmetric about the center 
of the galaxy will cause the particle to accelerate to a 
greater orbit for a given rotation velocity. Therefore, el- 
liptical orbits of particles in the Tk-r are inhibited and 
rtkr ~ 0. Remember one of the simplifying assumptions 
was the stars are stable in material composition rather 
than changing nis/m^. Changing msjrn^ would result in 
a slow inward spiral. Therefore, the radius rfrr of the 
orbits of the particles in the FRR become nearly circular 
and rfrr ~ 0. Thus, the mass in an elemental volume is 
Dfrr47rrfj.j.dr where Di„ is the density of matter at r^r. 
The Dfrr is approximately constant because the particles 
are of the same species. Therefore, the M in Eq. (|42() is 

Mfrr(rfrr) = ^Di„T^r\„ + MAfrr, (56) 

where 

MAfrr = Mtk -DfrrTrr^k (57) 

Mfrr(rfrr) is the mass within rfrr as a function of rtrr; 
Mtkrmax IS the total mass in the CR, Tc-k, KR, and 
Tk-r; and rtkrmax is the radius of the Tk-r. 

Substituting Eq. into Eq. and taking f 
yields 

Wfrr = ^-GDfrrTrrf^^ 
TOt 6 

m^rfrr 

-rfrr • ao(rfrr), (58) 

where Wfrr is the rotation velocity of a particle at rfrr. 

The r'^^ term was small at r = rkmax and is decreasing 
as rfrr increases. Therefore, the rj^ term is relatively 
much smaller than the term. Since the term for 
hydrogen differs from the rrig term for stellar material, 
the RCs will differ. By averaging the opposite sides to 
reduce the term as before yields 

TTh 4 

Wrr « — ^-GDfrrTrrl, -|-i^2, (59a) 

m, 3 

— -^frrslopc^frrmax ~^ -^frrintcrcept? (59b) 

where Xfrrsiopc and ii'frrintorccpt are the slope and inter- 
cept of the t'fj.j.-rfj.j. curve. 

Figured shows the Hi rotation curve of NGC 4321 0. 
It has the rotation velocity plotted well into the bulge. 
The curved line is in the FRR and is a plot of versus 
rfrr according to Eq. 
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FIG. 7: Plot of rotation velocity v^(10'' km^s~^) versus r 
(kpc) for NGC 4321 9] using the distance to the galaxy cal- 
culated using Cepheid variables to calculate r. A linear curve 
has been fitted to the SRR. A -r^ curve has been fitted to 
the FRR. 



2. Second RR 

The end of the FRR and the beginning of the SRR 
is caused by a change in the nis/m^ factor. The Space 
force starts to become significant because of ms/m^ rel- 
ative to the gravitational force. As the r increases, the 
particle orbits become rotationally flattened and remain 
circular pages 723-4]. The spherical symmetry is 
ended. If either the rotation velocity v^rr in the SRR, the 
Space force, or the gravitational force changes, the ac- 
celeration Tsn- in the SRR will change to restore balance 
and the circular orbits. Thus, the Principle of Negative 
Feedback applies and fsn- ~ 0. The orbits of the particles 
are nearly circular except for changes over a long period 
in the e; for changes in the elemental composition and 
mass of the particles such as by photon emission; and for 
changes in the r • term. Some of the galaxy's matter, 
especially the lighter matter with high nis/m^^ is moving 
radially. Most of a matter in the SRR is in a stable, radial 
position. Light material such as hydrogen and heavier el- 
ements differ in the mechanics defining their orbits due 
to the Space force. 

The mass in an elemental volume dr of the SRR can 
be modeled as a cylinder shell of height ffgrr -29.. page 
724] and with density Dsn- The mass in an elemental 
volume of the SRR is I?srr^^srr27rrsrrdrsi-r where Ts-c-c is 
the radius of a particle's orbit in the SRR. Unlike the 
inner regions, the rris/m^ of particles vary more slowly 
with increasing rgrr- The height profile of the bulge is 
a measure of the amount of mass versus m^jm^ particle 
type. The interplay of the terms of Eq. (|42|l determine 
the orbital radius of a particle. Thus, the mass Afgrr 



within the r^xx is 

Ms„ = (Arr-ffsrr27r) 4, + AfAr, (60) 

where 

A^Ar = Afft.„„ax - (L'srri?srr27r) r\„^^^; (61) 

''frrmax is the T to the outcr edge of the FRR; and Affrrmax 
is the mass in the CR, Tc_k, KR, Tk-r, and FRR. 

Inserting Eq. Ht)U|l into Eq. 142|) and with circular orbits 
yields 

-rsrr • Oo- (62) 

Since the particles in the SRR have larger ms/rrii, and 
since the r~l term is small at r = rkmax, the r~] term is 
small compared to the rsi-r term. Equation (|62|) suggest 
fg^j. is nearly independent of e and depends on galaxy 
parameters AfAn Dstr and i?srr- By averaging from side- 
to-side of the galaxy as before, the rgj-r • lo term is re- 
duced. Galaxies with other close galaxies in the v plane 
of rotation, close galaxies with high e values, or several 
close galaxies asymmetrically arranged, the residual of 
the side-to-side variation of the rs-cv • do term may be 
relatively significant. Asymmetry in the RC reflects this 
condition [l^. Therefore, 

Wg„ = i^srJ'srr + Ki^ , (63) 

where Kst and are the slope and intercept of the 
linear Ug^,. — rgrr relationship, respectively. The linear 
relationship of Eq. for NGC4321 is the straight fine 
plotted in Fig. |7| 

To test Eq. Ht)3|) , galaxies with a distance Dc calculated 
using Cepheid variables (the "Dz" of Freedman et al. 53]) 
and with a published Hi RC which included data points 
in the SRR were chosen. Figure |H1 shows the data points 
for each RC of the 15 selected galaxies. The straight lines 
in Fig.|Slare found by fitting a 0.99 or higher correlation 
coefficient, straight line to the first three data points im- 
mediately before the knee of the RC which indicates the 
transition region Tj-^s between the RR and SR. NGC 
4414 had only two data points before the knee. NGC 
3198 and NGC 0224 have several data points between 
the knee and the straight line of which no three are on a 
0.99 correlation straight line. 

The data for the 15 selected galaxies are shown in Ta- 
ble mil The morphology type data was obtained from 
the NED database. The morphology type code t and lu- 
minosity class code Ic data was obtained from the LEDA 
database [s^. The number "No." of data points es- 
tablishing the SRR, the values of the maximum radius 
Trmax (kpc) of the last data point on the SRR straight 
line, the maximum rotation velocity fimax (kms~^) of 
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FIG. 8: Plots of the rotation velocity v (10 km s ) versus radius r. The straight lines mark the SSR. 



the last data point on the SRR straight line, Kgr, and 
i^ir were obtained from the plots of Fig. |S1 The "Ref." 
lists the references for each RC. 

This sample has lovif surface brightness (LSB), medium 
surface brightness (MSB), and high surface brightness 
(HSB) galaxies; galaxies with a range of the 21-cm line 
width W20 at 20 percent of peak of from 120 km s~^ 
to 607 kms~^; includes LINER, Sy, HII and less active 
galaxies; galaxies which have excellent and poor agree- 
ment between distances Dtf calculated using the TuUy- 



Fisher relation (TF) and Dc, a Dc range of from 0.75 
Mpc to 17.39 Mpc; field and cluster galaxies; and galax- 
ies with rising, flat, and declining RCs in the SR. 

The CUM predicts a vfj.,. — r^.^. relationship if a galaxy 
has spherical symmetry in a FRR. At least five of the 
galaxies in Fig. |S1 have data points sufficiently far inward 
to show the vf^^ — rf^^ relation and don't. The CUM 
suggests these galaxies lack a FRR (spherical symmetry) . 

As Tsrr iucrcascs, rris/m^ increases and the predomi- 
nant particles have less mass and decreasing density. De- 
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TABLE III: Data for the chosen galaxies. 



„ . 


ixiorpnoio^iy 


t' 




RC^ 


D 


No. 














Kir'' 


Rcf. 


"M" f~i m o /I 


Ort-I^S ) D 


3 


2 


r 


0. 79 


3 


10.3ib3 


1 




- 

5 


lOOit 


200 


- 

-5 


900±11 


000 




"\Tf~T~' mnn 
IN y~-r\^ UoUU 




6.9 


6 


±1 


2.00 


6 


3. Iit2 


2 


72ib 13 


1 


700± 


100 




400± 


400 




IN i-jV^ uoyo 


Ol\-\ S / CQ 


D 











2 4ih 1 








400± 


300 




000± 4 


000 






Or\.D\S J CQ 


D 





p 


3 22 


o 
O 




Q 


124ibl0 


2 


860± 


80 


2 


000± 1 


000 






bA( rJb;ijiiNrL/H by 


3 


1 


F 


14.07^ 


3 


8.3ib2 


4 


324±18 


26 


400± 


700 


-8 


000± 4 


000 


[57. 581 




9 A f^'j'lah-T JNF.R R 


2.4 


2 


j3 


3 63 




7.4±0 


g 


243il4 


3 


900± 


200 


30 


000± 5 


000 




NGC 3198 


SB(rs)c 


5.2 


3 


F 


13.80 


4 


5. Oil 


2 


134±15 


3 


100± 


300 


2 


000± 3 


000 




NGC 3319 


SB(rs)cd;HII 


5.9 


3.3 


R 


13.30 


6 


8.7±1 


7 


101± 5 


1 


300± 


100 


-1 


000± 4 


000 




NGC 4258 


SAB(s)bc;LINER Syl.9 


4 


3.9 


D 


7.98 


3 


2.3±0 


7 


197±26 


36 


000±10 


000 


-43 


000±30 


000 




NGC 4321 


SAB(s)bc;LINER HII 


4 


1 


D 


15.21 


6 


5.8±0 


8 


233±14 


5 


100± 


500 


20 


OOOilO 


000 




NGC 4414 


SA(rs)c? LINER 


5.1 


3.6 


D 


17.70 


2 


6.4±2 


9 


247±45 


6 


000 




20 


000 




[63 




NGC 4535 


SAB(s)c 


5 


1.9 


D 


15.78 


3 


4.6±1 


2 


166±12 


7 


400± 


300 


-7 


000±17 


000 






NGC 4548 


SBb(rs);LINER Sy 


3.1 


2 


D 


16.22 


3 


10.6±2 


9 


195±12 


3 


700± 


300 


-1 


000±17 


000 


165 




NGC 5457 


SAB(rs)cd 


5.9 


1 


D 


6.70 


4 


8.5±1 


3 


228±12 


4 


000± 


400 


20 


000±20 


000 


m 




NGC 7331 


SA(s)b:LINER 


3.9 


2 


F 


14.72 


3 


2.0±0 


6 


211±20 


11 


000± 2 


000 


23 


000±14 


000 


w 





"Galaxy morphological from the NED database. 

''Galaxy morphological type code from the LEDA database. 

'^Galaxy luminosity class code from the LEDA database. 

■^Galaxy's HI rotation curve type according to slope in the outer 
SR region. R is rising, F is fiat, and D is decl inin g. 

•^The distance Dc (Mpc) to the galaxy from l53l unless otherwise 
noted. 

■'^The maximum rotation velocity Urmax (kms— 1), the radius 
''rmax (kpc) of i;rmax, and the No. of points on the straight line 
from the curves in Fig. |H] The error of rrmax and Vsn is ±5% of 
the value plus the difference of the chosen data values and the next 
data point value. 

f Least squares slope of - for each galaxy in the SRR of HI 
rotation curve (km^ s~^ kpc"'^). 

''Least squares intercept of - for each galaxy in the SRR of 
HI rotation curve (km^ s~^). 

'The Dc is from 6^. 



creasing amount of mass in a given cylindrical shell and 
the rotation flattening effect causes iJgrr to decrease to 
a disk. This is end of the RR and the beginning of the 
transition region Ti._s between the RR and SR. 

3. RR data 

Because this change is caused by the change in elemen- 
tal types which depends on the e of a galaxy, the end of 
the RR is occurring under similar conditions in all galax- 
ies. Thus, the mass Mj-max of the galaxy inside a sphere 
with a radius of rrmax can be used to compare parameters 
among galaxies. 

Evaluating Eq. H42|l at rrmax and rearranging terms 
yields 

rrmaxt^rmax = — GA/rmax -GsK^dr^d + Km, (64) 

where 

Kii-Y ^rmax^i'niax rrmax^rmax * ^o- (^^) 

The rrmax^rrnax term is proportional to the effective 
mass McSr inside a sphere with a radius of rrmax- Since 

' rmax 

~ 0, Kirr is a relatively small value. However, in 
addition to the close-galaxy caveat, a very young galaxy 
or a galaxy with significant infall of matter such as col- 
liding galaxies can cause a deviation in Kij-i-. 



Since rrmax is in the RR and not in Tr_s or SR, Wrmax 
can be found in rising, flat and falling RCs. The gas 
spectral lines such as Hq reflects the motion of a hot 
star and not hydrogen gas. The elemental content of 
hot stars changes with r which changes the Ws/mi term. 
Therefore, "yrmax changes with each star elemental type. 
The Hi RC is monitoring the same particle type in each 
galaxy from the RR to the SR. Therefore, the Hi RC was 
used to test Eq. (|^ . 

The red data listed in Table HVl were calculated from 
Dc in Table ITTll and from the average angular (observed) 
radius ra (arcsec). The ra for each sample galaxy was 
measured from the Hubble Space Telescope (HST) im- 
ages listed in Table HVl The Visual Basic program listed 
in Appendix ^ was used to extract the I ~ r curve and 
calculate ra. The selection of each image was based upon 
being in the V band, upon having a sufficient area of the 
center of the galaxy in the image so ra on both sides 
of the galaxy could be measured, and with sufficient ex- 
posure time to have the I value at ra sufficiently large. 
The I — r curve obtained from the HST images using the 
WFPC2 instrument are shown in Fig. M NGC 2841 had 
only one side of the galaxy in the image. 

The straight line in each I — r curve of Fig.|51is drawn 
between the pixels the Visual Basic program selected to 
calculate ra. Note, the line is at differing / levels among 
the galaxies and the slope of the line indicates asymmetry 
in / value for the calculation of ra in the same image. 
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TABLE IV: Data of the galaxies in Fig. ^ 



Galaxy" 


HST 






line'^ 




image 


filter 


pc 


n m p s 


N0224 


U2LG0209T 


F555W 


2.87±0.38 


- 6 7 7 


N0300 


U6713701M 


F547M 


3.44±0.97 


- 2 5 3 


N0598 


U3MR0106M 


F555W 


1.50±0.41 


- 3 6 5 


N2403 


U6712504R 


F547M 


6.61±1.56 


- 3 4 4 


N2841 


U65M7008M 


F555W 


15.68±6.82 


- 5 4 5 


N3031 


U32L0105T 


F547M 


5.11±1.76'' 


3 5 5 6 


N3198 


U29R1602T 


F606W 


64.20±6.69 


- 1 1 1 


N3319 


U29R1702T 


F606W 


19.05±6.45 


- 2 4 2 


N4258 


U6712101R 


F547M 


13.41±3.87'' 


2 3 2 4 


N4321 


U2465107T 


F555W 


23.66±7.38 


- 3 3 4 


N4414 


U67N8102R 


F606W 


39.95±8.58 


- 3 2 3 


N4535 


U2782603T 


F555W 


31.25±7.65 


- 2 2 3 


N4548 


U34L0303T 


F555W 


32.34±7.87 


- 3 3 3 


N5457 


U6712704R 


F547M 


5.32±3.25 


- 5 5 6 


N7331 


U41V1802M 


F450W 


19.23±7.14 


- 2 1 4 



""I" means IC and "N" means NGC. 

''The rj(j was measured as outlined in tlie text. Tlie error is ± 0.1 
arcsec. 

'^Tlie line is tlie "n" listed in Table|n]and shown in Fig.|S] the "m" 
listed in Table Ivl and shown in Fig. 1101 the "p" listed in Table IvTI 
and shown in Fig. 1111 and the 's" listed in Table fvTTl and shown in 

Fig. EH 

''The values differ from Table because of differing distances 
used. 



TABLE V: Data of the lines in Figures [TUl 

line'' Corr.<^ F Test Kstt"^ f Ki„'= f 

lO'^km^s-^ IQ-lkpckm^s-^ 

1 filled diamonds ND^ ND OTi 

2 filled squares 0.99+ 0.99+ 0.42+0.03 0.3+0.5 

3 filled triangles 0.97 0.94 1.1 +0.1 -0.4+3 

4 ND 

5 filled circles 0.97 0.96 5.2 +0.9 9 +8 

6 open diamond ND ND 16.8 



"The integer denoting the place of the line in the order of increas- 
ing Xsoft. 

''The identification of the line symbol in Figurcs ^lOl 
'^The correlation coefficient. 

''The least squares fit slope of the ''rmaxi'rniax " '"ed lines. 
"^The least squares fit intercept of the rrmaxt^rmax " ^ed lines. 
•''The lines are calculated using (r^d 

'^rmax) = (0:0) ^ ^ '^^'t^ 

point. 

f "ND" means no data because there are only two points (one data 
point). 



Therefore, the selection of ra is not based on selecting an 
isophotal value. The ra is not an isophotal based effective 
radius. 

Figure ITUl is a plot of rj-maxWrmax versus r^d- Because 
the sample is small for six galaxy line classifications (four 
of the six line classifications have less than four galax- 
ies), the (r-edj'^rmaxi'rmax) = i^fi) Point was included in 
all classifications. Table Ivl lists the calculated values for 
the hues of Fig. Uni 

The lines are straight lines. Therefore, the equation 



for each line is 

^rmax^rmax — -^sir^ed -^irr; (66) 

where ifgrr and Xirr are the slopes and intercepts of the 
lines in Fig. ^1 respectively. 

Since Mrmax > 0, the intercepts of the lines in Fig. [TUI 
nearly intersect at (red,rrmaxi'?max) = (0>0); and Ki„ is 
determined by neighboring galaxies, 

Mrmax = K^-^C = K^^K^dTed, (67) 

where K^^ is the proportionality constant. 

The slopes among the lines in TableElobey the relation 

l°glO ( , ^2" o ) = -^sni-m + A'inr , (68) 

Vkm s^^/ 

where K^^, = 0.40±0.02 and K^^, = 5.79±0.06 at 1 cr are 
the least squares fit of the slope and intercept of the linear 
relation of Eq. (|68|) . respectively. The m is an integer 
of one to six depending on the line from Table The 
correlation coefficient and F test for Eq. H68|l are 0.99 and 
0.99, respectively. Note, the lowest slope corresponds to 
m=I in Eq. (|68|l and the lines nearly intersect at a point. 
The KsvT is nearly independent of the distance used to 
calculate the rrmaxf^max a.nd r^d if the same distance is 
used for both. The values of logio(e)Ri 0.434 is within 
1.7 cr of iCsm- 

For the rnnaxWrmax ~ ''ed relationship, 

TTl TTi 

Ksrr = ^GK,,K,d - — Gsif.d. (69) 

TO,, 

Therefore, restating Eq. H68|l as a strong Principle of 
Repetition yields 

Ksn^i^scre", (70) 

where K^or = (4.8 ± 0.9) x 10^ km^ s"^ at I a. 
Combining Eqs. H66|) and H70|l yields 

'"rmax^rmax — -^scr ^ ^ed; ('^-i) 

where Ki„ = 0. 

The standard deviation of the function are the error 
bars of r^d hsted in Table |IV| The = H and w 
70% between Eq. H71|l and the measured rrinaxi'rmax (^ot 
including the (r<;d, rrmaxi^?max) = (O^O) point). Since 
is greater than the 5% significance level, this test does 
not invalidate Eq. H71|l . The ''rmax^^rmax i^ a function of 
r^d and to. The r^d is a parameter of each galaxy, the 
TO is a galaxy classification parameter, and the e implies 
the Principle of Repetition is applicable to the ri-max^'rmax 
parameter, which is proportional to the effective mass at 

^rmax ■ 

Since rrmax^'rmax i^ ^t the outer edge of a TOs/to^ strata, 
it may be used to compare parameters among galaxies. 
Posit that Ti-max and w^max are related to r^d, also. 
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FIG. 9: Plots of the I — r profile for the sample galaxies. The straight line is drawn between the chosen pixels. The slope on 
the line indicates the asymmetry. The horizontal scale for each plot is approximately ten times the r^. 



Figure ITTI shows a plot of rj-max versus r^d for the sam- 
ple galaxies in Tables IIIII and IIVI Because the sample values are shown in Table IVII 
is small for seven galaxy integer p classifications (all of 
the seven p classifications have less than four galaxies), 
the (r^d, ''rmax) — (O7O) point was included in all classifi- 
cations. The linear relationships are 



Ted relationship, respectively, and the respective 
ire shown in Table IVTI 
The relationship of the -fiTsrd values is 

logio -^sid = KscidP + Kicrd, (73) 



(72) 



where K^crd = 0.28 ± .02 and Xicrd = 1-63 ± 0.08 at 1 a. 
Note, logio2 « 0.3010. The p is a power of two rather 
where Ksrd and Xird are the slope and intercept of the than e. This implies the strong Principle of Change is 
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FIG. 10: Plot ol the measured Lax (lO^kpC Ws"^) 

versus the distance red (kpc) along the major axis of a HST 
image to the first discontinuity or inflection point using the Dc 
stated in Table HV1 The filled diamonds, filled squares, filled 
triangles, filled circles, and open diamond are data points 
along a linear rrmax^rmax - '"ed line using data from Tables UTTI 
and IIVI and correspond to m = 1, m = 2, m — S, m = 5, 
and m = 6 as listed in Table IVI respectively. The m=6 line is 
drawn between NGC 0224 and (rEd,»'rmaxi'?max) = (0,0). The 
error bars were omitted for clarity. 
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FIG. 11: Plot of the measured rrmax (kpc) versus the dis- 
tance Tjd (kpc) along the major axis of a HST image to the 
first discontinuity or infiection point using the Dc stated in 
Table EVl The filled diamonds, filled squares, filled triangles, 
filled circles, open diamonds, open square, and open triangle 
are data points along a linear rrmax - r^d line using data from 
Tables UTTI and Hvl and correspond top=l,p = 2,p = 3, 
p = 4, p = 5, p = 6, and p = 7 as listed in Table IvTI respec- 
tively. The lines are drawn to include the point (r^d, rrmax) = 
(0,0). The error bars were omitted for clarity. 





TART P VT- 


Data of the lines 


ill i^Uicb 


rm 

[ill 


pa 


line' 


Corr."^ 


F Test 
















kpc 


1 


filled diamonds 


0.99+ 


0.99+ 


76 ± 9 


0.2 ±0.4 


2 


filled squares 


0.99+ 


0.99+ 


156 ± 8 


0.03±0.2 


3 


filled triangles 


0.98 


0.98 


310 ± 60 


-0.3 ±1 


4 


filled circles 


0.98 


0.98 


459 ±60 


0.6 ±0.8 


5 


open diamonds 


0.96 


0.95 


1500 ±300 


-0.6 ±1 


6 


open square 


ND9 


ND 


1620 





7 


open triangle 


ND 


ND 


3590 






"The integer denoting the place of the line in the order of increas- 
ing A'srd- 

''The identification of the line symbol in Figures ^lll 
"^The correlation coefficient. 

''The least squares fit slope of the rrmax - r^d lines. 
■^The least squares fit intercept of the rrmax - r^d lines. 
^The lines are calculated using (r^d, rrmax) = (0,0) as one of the 
data points. 

9 "ND" means no data because there are only two points (one data 
point). 



applicable and 



srd 



scrdt 



2^ 



(74) 



where K^, 



rdt 



36 ± 8 at 1 cr. 



Combining Eqs. 
yields 



((7^ and ((73Il and setting K„ 



scrdt 



(75) 



The standard deviation of the \^ function are the error 

10 and Pv ~ 76% 



bars of r^d hsted in Table llVl The x 
between Eq. H75|l and the measured Trmax (not including 
the (redj^'rmax) — (0,0) point). Since is greater than 
the 5% significance level, this test does not invalidate 
Eq. (|75|1 . The rj-max is a function of r^d and p. The r^d is 
a parameter of each galaxy, the p is a galaxy classification 
parameter, and the 2 implies the Principle of Change is 
applicable to the rrmax parameter. 
Posit 



irvd ; 



(76) 



where -fiTsrvd and i^n-vd are the slope and intercept of the 
linear relation, respectively. 

Using the data from Tables HTTl and HVI the plot of Wrmax 
versus r^d is shown in Fig. 1121 and the data for the lines is 
listed in Table IVIll Because the sample is small for seven 
galaxy integer s classifications (five of the seven s classi- 
fications have less than four galaxies), the (j^ediWrmax) — 
(0,0) point was included in all classifications. The rela- 
tionship of the slopes of the lines is 



log 



10 



-^srvd 



km^ s~^kpc ^ 



scrvd 



icrvdi 



(77) 



where ivT^crvd = 0.31 ± .02 and iCicrvd = 5.16 ± 0.07 at 1 
a. 

The strong Principle of Change is applicable and 



■vd 



scrds 5 



(78) 
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FIG. 12: Plot ol the measured Urmax (km'^s"^) versus the 
distance red (kpc) along the major axis of a HST image to the 
first discontinuity or inflection point using the Dc stated in 
Table Eni The flUed diamonds, filled squares, fiUed triangles, 
filled circles, open diamond, open triangles, and open square 
are data points along a linear w^max - '"ed line using data from 
Tables ITTTl and HVI and correspond to s = 1, s = 2, s = 3, s = 4, 
s = 5, s = 6, and s = 7 as listed in Table IVITl respectively. 
The error bars were omitted for clarity. 



TABLE VII: Data of the lines in Figures [HI 







Corr."^ 


F Test 


-^srvd ^ 

10^ kpckm^ s-2 


10^ km^s-^ 


1 


filled diamonds 


NDS 


ND 


0.28 





2 


filled squares 


ND 


ND 


0.54 





3 


filled triangles 


0.95 


0.92 


1.3 ±0.3 


-1 ± 7 


4 


filled circles 


0.99 


0.98 


2.3 ±0.2 


1 ± 3 


5 


open diamonds 


0.99+ 


0.99+ 


6.8 ±0.1 


-1 ± 1 


6 


open triangles 


0.99 


0.97 


10.6 ±0.04 


0.2± 4 


7 


open square 


ND 


ND 


16 






"The integer denoting the place of the line in the order of increas- 
ing Ksrvd- 

'The identification of the line symbol in Figures |l2l 
'^The correlation coefficient. 

"^The least squares fit slope of the f^nj^x " ''ecI lines. 
'^The least squares fit intercept of the v^^^-^ - r^^ lines. 
.'^The lines are calculated using {rf:^,v'^^^^) = (0,0) as one of the 
data points. 

9 "ND" means no data because there are only two points (one data 
point). 

where Kgervds = (1.6 ±0.3) x 10^ km^s~^kpc~^ at 1 a. 
Combining Eqs. and yields 

^rmax — -^scrvds 2 T'gcj, ('^'^) 

where -fi'irvd = 0. 

The standard deviation of the function are the er- 
ror bars of r^d hsted in Table E| The = 8.2 and 
~ 88% between Eq. (|7^ and the measured v^^^^ 
(not including the (f'ed, i^rmax) = (0,0) point). Since 



is greater than the 5% significance level, this test does 
not invahdate Eq. l(7^ . The w^max ^ function of r^d 
and s. The r^d is a parameter of each galaxy, the s is 
a galaxy classification parameter, and the power of two 
implies the Principle of Change is applicable to the w^max 
parameter. 

A feature of Figs. E| and ^1 is the lack of data 
points in the upper right of the plots (galaxies with rel- 
atively high Ted and high r„nax, W?max, and r„„axWrmax 

values). 

E. SR 

In the bulge, strong gravitational and Space forces sep- 
arate the regions into strata by particle ms/m^ type. In 
the Ti._s, the decreasing height i?trs of a galaxy's cylin- 
drical shells results in a lower slope of the plot of the 
square of the rotation velocity v^^.^ in the Ti._s versus 
radius rtrs in the Tj—s (^^tis " ''trs curve) than in the RR. 

In the SR, the AI and e forces are less than in the KR. 
Therefore, the tendency to separate into separate ms/m^ 
regions is less than in the KR. Therefore, the regions 
are less distinct. The SR could be modeled with these 
independent regions as rings as is often done to analyze 
rotation curves. Since the change of ms/m^ from one ring 
to the next is small (unlike in the bulge), the form of the 
fundamental equation of motion for each ring will change 
little. Also, since the elemental composition of the stars 
is changing, the border of the rings may be blurred. For 
purposes herein, to calculate the mass Ms within a sphere 
with a radius rg in the SR, the SR is considered as a whole 
with a smoothly varying rather than stepwise changing 
ms/m,. 

The thin disk of the SR is modeled as a cylinder of 
density Dg and height Hg, then the Ms is 

Ms = {DsHs2n) Ts + Mas, (80) 

where 

Mas = Mtrs ~ {DsHs2Tr) r^max (81) 

and Mtrs is the mass inward of the SR. 

Because the rotation appears as a rigid body, the pro- 
jection of the Ms into the RR may be real and not just 
a mathematical convenience. 

Inserting Eq. H80|l into Eq. (|42|) and having circular 
orbits [2^ page 725] yields the rotation velocity Vs in the 
SR 

«2 = (j^GDsHs27r^ Ts 

_^ m,GM^s-n.sGse _^^^^^_ (82) 
m,rs 

As in the RR, the term is nearly constant. By 
averaging from side-to-side of the galaxy as before, the 
rs • flo term is reduced. However, unlike in the RR, 
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DgHs < DsirHsrr. Therefore, rg • is more signifi- 
cant. In cases where there are other close galaxies, close 
galaxies with moderate e values, several close galaxies 
asymmetrically arranged, or a close galaxy with a high 
e, the residual of the side-to-side variation in the • Qq 
term may be relatively significant (see Hodge and Caste- 
laz Therefore, 



where 



Kssrs + Ks 



(83) 



where Kss and i^^si are the slope and intercept of the 
linear relationship, respectively. 

By reasoning similar to the SRR region, the maximum 
rotation velocity Wsmax in the SR is also a linear function 
of e. Therefore, for the intrinsic galaxy with rg • flo ~ 



where 



seslopc 



(84) 



(85) 



and the constant K^s — -^^smax/''ed- 

However, unlike the RR, the last rotation velocity data 
point measured may not be Wsmax and will be strongly 
influenced by • Oq. 

If the impact of the rg • Oo is small, the RC in the 
SR is linear and rising. Another result of a small Tg • 
Oo is that the difference of rotation velocity at a given 
radius from one side to the other of a galaxy is small. 
However, if rg • has a value comparable to the other 
terms of Eq. (|83|) and V^(rs»ao) ~ 0, the rs'flo term also 
increases approximately linearly. Therefore, the SR RC 
will be approximately flat with a slope depending on the 
relative value of the terms. If V^(rs • flo) is significant, 
the average rotation curve will decline and there will be 
significant asymmetry in the RC at a given rg. Since 

• flo is in the plane of the galaxy, the asymmetry due 
to neighboring galaxies found by Hodge and Castelaz ^| 
is in accordance with Eq. H83|) . 

Consider two adjacent ring volumes in the disk with 
different elemental composition between the rings and 
with uniform elemental composition in each ring. Within 
each ring, Eq. H44|l applies and the RC will be Keplerian. 
Because of the differing elemental composition, the fg 
between the rings obey Eq. H83(l and the outer ring will 
have a greater Ug than the inner ring if rg • — 0. 
However, the differing ms/m^ is small. Therefore, the 
detection of the slight decline within each ring may be 
difficult. 



DR 



The DR is where particles are no longer rotating 
around the galaxy. The galaxy forces are comparable 
to the forces from other galaxies. Equation (|42|) reduces 
to 



Clod 



all galaxies 

E 

Sinks 



mgGAIj — rUsGsej 



fe=l 



(87) 



and the sum is now over all galaxies. 

Equation H87() has the form of the Newtonian gravi- 
tational force with a anti-gravitational like force, both 
centered at the center of a galaxy. 



IV. OTHER RELATIONSHIPS 

A. SAMPLE DATA 

The sample of gala xies was drawn from the sample of 
Freedman et al. |53l |. who present a sample of 31 spi- 
ral galaxies with distances Dc (Mpc) calculated using 
Cepheid variables ("i'z" in Freedman et al. (53l|). Since 
Ted was desired, the sample required each galaxy to have 
a HST, optical FITS image of the central few kpc of 
the galaxy. NGC 0925 lacked an adequate HST image 
and, therefore, was excluded from the sample. NGC 2841 
'68'! was added to the sample since the Dc determining 
method was consistent with Freedman et al. |53| and be- 
cause NGC 2841 is a potential falsifier of the MOND 
model 5]. 

The data for the 31 selected galaxies are shown in Ta- 
bles Ull 



rdrd - fd • ctod, 



(86) 



Villi and IIXI The morphology type data were 

obtained from the NED database. The luminosity class 
code Ic and de Vaucouleurs radius R25 data were obtained 
from the LEDA database. The Hi (21 cm) line width W20 
(kms~^) at 20% of the peak value corrected for inclina- 
tion data were calculated from "incl" and "w20" param- 
eters from the LEDA database. The rrih = "a magni- 
tude in the visible range" - "extinction" from the NED 
database. The methods and equations used to calcu- 
late the Fp, s, u, w, and x values are described in Sec- 
tion |IVB1 

The red data listed in Table Ivml were calculated from 
Dc listed in Tables UTTknd IVnTl and from r^. The for 
each sample galaxy was measured from the HST images, 
which were obtained using the WFPC2 instrument, listed 
in Table IVlTll The Visual Basic program listed in Ap- 
pendix B was used to extract the I—r curve and calculate 
ra. The selection of each image was based upon being in 
the V band, upon having a sufficient area of the center 
of the galaxy in the image so Ta could be measured, and 
upon the / value at being sufficiently large. Figure |3 
shows the I — r curves for the galaxies with an RC. The 
remaining I — r curves are shown in Fig. 1131 The straight 
line in each I — r curve of Fig. 1131 is drawn between the 
pixels the Visual Basic program selected to calculate r^. 
Note, the line is at differing / levels among the galaxies. 
The slope of the line indicates asymmetry in / value for 
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TABLE VIII: Data for the chosen galaxies without rotation curves. 



Galaxy 


morphology" 


hst' 














image 


filter 




Mpc 


xlO"^kpc 


11613 


IB(s)m 


U4040201R 


F555W 


9.8 


0.65 


1.2± 0.3 


14182 


SA(s)m 


U6G24404R 


F606W 


8.1 


4.49 


5.0± 2.2 


N1326A 


SB(s)m: 


U34L3401O 


F555W 


7.8 


16.14 


14.4± 7.8 


N1365 


(R')SBb(s)b Syl.8 


U2KV010XT 


F547M 


1 


17.95 


30.5± 8.7 


N1425 


SA(rs)b 


U29R5P02T 


F606W 


3.1 


21.88 


31.9±10.6 


N2090 


SA:(rs)b 


U29R0W02T 


F606W 


3.9 


11.75 


19.8± 5.7 


N2541 


SA(s)cd LINER 


U29R0Z02T 


F606W 


6 


11.22 


9.1± 5.4 


N3351 


SB(r)b;HII Sbrst 


U67N3002R 


F555W 


3 


10.00 


17.4± 4.9 


N3368 


SAB(rs)ab;Sy LINER 


U6EAN504M 


F606W 


3 


10.52 


22. 2± 5.1 


N3621 


SA(s)d 


U29R1I02T 


F606W 


5.8 


6.64 


10.0± 3.2 


N3627 


SAB(s)b;LINER Sy2 


U29R1K02T 


F606W 


3 


10.05 


55.3± 4.9 


N4496A 


SB(rs)m 


U2690102T 


F555W 


5.9 


14.86 


14.3± 7.2 


N4536 


SAB(rs)bc HII 


U2DT0601 


F555W 


1.9 


14.93 


114.1± 7.2 


N4639 


SAB(rs)bc Syl.8 


U2NU0202P 


F555W 


3.8 


21.98 


135.2±10.7 


N4725 


SAB(r)ab pec Sy2 


U67N4602R 


F606W 


2 


12.36 


28. 1± 6.0 


N5253 


Im pec;HII Sbrst 


U2E68101T 


F606W 


ND 


3.15 


3.9± 1.5 



"Galaxy morphological from the NED database. 
'Galaxy Hubble Space Telescope image from MAST. 

'^Galaxy data from the LEDA database. 

"^The distance Dc (Mpc) to the galaxy from Freedman et al. l53l 
unless otherwise noted. 



TABLE IX: Additional data for the chosen galaxies. 



Galaxy 






^25° 


F ' 




S U W X 




km s ^ ^ 




kpc 


xio" 


mag. 




Galaxies with rotation curves 














N0224 


548±10 


21 


.4±1.4 


ND' 


4.26 


7 8-9 


N0300 


233±11 


5 


.7±1.6 


3.22±0.93 


8.95 


3 6 5 5 


N0598 


243 ± 8 


8 


.1±0.5 


ND 


6.09 


5 8-8 


N2403 


279±18 


11 


.Oil. 8 


2.98±0.86 


8.77 


4 6 5 6 


N2841 


655± 7 


15 


.5±2.0 


1.04±0.30 


10.09 


5 5 3 7 


N3031 


517± 4 


11 


.8±1.0 


0.98±0.28 


7.75 


6 6 4 8 


N3198 


343 ± 7 


15 


.6±2.6 


2.89±0.84 


10.87 


13 2 4 


N3319 


263±13 


11 


.4±1.9 


2.01±0.61 


11.48 


2 4 3 5 


N4258 


461± 5 


19 


.7±3.3 


3.68±1.07 


9.10 


4 6 4 7 


N4321 


540±10 


16 


.8±2.5 


14.90±4.32 


10.00 


4 5 5 7 


N4414 


509 ± 9 


9 


.1±2.0 


1.72±0.50 


10.94 


3 3 2 5 


N4535 


423±12 


15 


.9±1.4 


6.41±1.86 


10.59 


3 4 4 6 


N4548 


449±31 


12 


.4±1.1 


0.81±0.23 


10.90 


3 4 2 5 


N5457 


502±12 


30 


.1±6.7 


3.14±0.91 


8.31 


6 7 5 9 


N7331 


538± 7 


22 


.9±3.4 


1.00±0.29 


10.02 


4 5 3 7 


Galaxies without rotation curves 














11613 


68± 6 


1 


.5±0.3 


ND 


9.86 


15-2 


14182 


134±11 


3 


.9±0.3 


22.80±6.61 


13.00 


15 7 1 


N1326A 


130±16 


4 


.0±0.9 


24.90±7.22 


13.77 


13 6 3 


N1365 


465±15 


28 


.6±2.5 


1.66±0.48 


10.32 


3 5 3 7 


N1425 


392±10 


19 


.2±2.5 


3.48±1.01 


11.87 


3 4 3 5 


N2090 


319±10 


9 


.8±2.9 


1.60±0.46 


11.99 


2 4 3 4 


N2541 


230±11 


9 


.4±2.4 


4.19±1.22 


12.09 


2 5 5 5 


N3351 


424± 6 


11 


.0±0.7 


1.55±0.45 


10.49 


4 4 3 5 


N3368 


439±14 


11 


.9±1.0 


32.70±9.48 


10.06 


4 4 6 6 


N3621 


306± 6 


11 


.3±2.5 


0.05±0.01 


9.79 


3 5 16 


N3627 


449±10 


12 


.7±1.4 


0.70±0.20 


9.64 


2 3 15 


N4496A 


235± 7 


8 


.2±0.5 


33.20±9.63 


11.93 


2 4 6 5 


N4536 


394±14 


15 


.7±1.7 


19.50±5.66 


11.16 


12 4 3 


N4639 


420±43 


9 


.7±0.6 


26.90±7.80 


12.19 


114 2 


N4725 


510±14 


19 


.7±1.3 


19.70±5.71 


10.08 


4 5 5 6 


N5253 


101±14 


2 


.2±0.2 


0.58±0.173 


10.68 


14 4 4 



"Galaxy data from the LEDA database. 
'The units of Fp are ergcm"'^ Mpc~^. 

'^"ND" means No Data because the Milky Way would be one of 
the close galaxies. 
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the calculation of ra. Therefore, the selection of ra is not 
based on selecting an isophotal value. The ra is not an 
isophotal based effective radius. 

This sample has low surface brightness (LSB), medium 
surface brightness (MSB), and high surface brightness 
(HSB) galaxies; galaxies with a range of W20 from 68 
km s~^ to 655 kms~^; includes LINER, Sy, HII and less 
active galaxies; a Dc range of from 0.65 Mpc to 21.98 
Mpc; field and cluster galaxies; and galaxies with rising, 
flat, and declining RCs in the disk region. 

B. RESULTS 

Figure ini shows a plot of W20 versus Urmax- The filled 
diamonds use the measured values of W20 and frmax of 
the galaxies with an RC. The line is a plot of. 




where Ki = 7.3 and K2 = 0.78. 

Equation (|79|l relates r^d and Wrmax- The s for each 
galaxy in Table Hxl without a RC was chosen to provide 
the data point closest to the line of Eq. (|88|l . 

The uncertainties of r^d listed in Table were used to 
calculate the standard deviation of the function. The 
= 8 and w 86% between Eq. (|S5|l and the mea- 
sured W20- Since is greater than the 5% significance 
level, this test does not invalidate Eq. The W20 is 

a function of Wrmax that is a function of r^d and s. The 
Tjd is a parameter of each galaxy, the s is a galaxy clas- 
sification parameter, and the "2" implies the Principle of 
Change is applicable to the W20 parameter. 

The W20 is used by the Tully-Fisher relationship |6^ 
ItoI I to calculate total absolute magnitude (mag.) in 
band /3 

=i^4/3+if5/3 log(W^20), (89) 

where K^p and K^p are constants that depend on the 
band under consideration. Therefore, is a function 
of s and r^d- 

The correlation coefficients of W20 versus rrmax and 
''rmaxWrniaxJ which is proportional to the effective mass 
within a sphere of radius rj-max, are 0.38 and 0.80, re- 
spectively. 

Figure [TBI is a plot of R25 versus r^d- Each line in 
Fig. El is a different "classification" . Because four of 
the eight classifications had less than four galaxies, the 
(?'£d,^25) = (0,0) point was included in all classifications. 
Table IXI lists the calculated values for the lines of Fig. 1151 

The lines of Fig. E| are straight lines. Therefore, the 
equation for each line is 

^25 = Ksr25red + Ku2b, (90) 

where ii'sr25 and K]j-25 are the slopes and intercepts of 
the lines in Fig. ^1 respectively. 



TABLE X: Data of the lines in Figures [T^ 





Ime" 


Corr."^ 


F Test 




Ki 




1 


filled diamond 


NDf 


ND 


71 







2 


filled square 


ND 


ND 


138 







3 


filled triangles 


0.99+ 


0.99+ 


233± 9 


1 


+ 1 


4 


filled circles 


0.95 


0.87 


490± 60 


1 


+ 1 


5 


open diamonds 


0.95 


0.87 


800+100 


1 


+2 


6 


open squares 


0.97 


0.95 


1400+200 


1 


+2 


7 


open triangles 


0.99+ 


0.99+ 


5680+ 80 


0.2+0.3 


8 


open circle 


ND 


ND 


7450 








"The integer denoting the place of the line in the order of increas- 
ing J^sr25- 

'The identification of the line symbol in Figures ^lSI 
'^The correlation coefficient. 

''The least squares fit slope of the iJ25 - ^ed lines. 
■^The least squares fit intercept of the R25 - r^^ lines. 
^The lines are calculated using {r^i,R2h) = (0,0) as a data point. 
9 "ND" means no data because there are only two points (one data 
point). 

Since Kn:2T, ~ 0, the slopes among the lines in Table Ixl 
obey the relation 

logio (^^) =i^suU + i^iu, (91) 

where K^^ = 0.30 ±0.02 and i^iu = 1.52 ±0.08 at Icr are 
the least squares fit of the slope and intercept of the linear 
relation of Eq. H91|l . respectively. The u is an integer of 
one to eight depending on the line from Table H The 
correlation coefficient and F test for Eq. H91|l are 0.99 and 
0.98, respectively. Note, the lowest slope corresponds to 
u = 1 in Eq. 1191(1 . The i?25 1 Ted is nearly independent of 
the distance used to calculate the i?25 and r^d if the same 
distance is used for both. The value of login (2)« 0.301 
is well within la of i^su- 

Therefore, restating Eq. 1(9111 as a strong Principle of 
Change and combining Eqs. (I90II and ((91() yields 

i?25 = ifRr2"r,d, (92) 

where — (32 ± 7) at la and K„25 = 0. Note the lack 
of data points in the upper right of Fig. 1151 

The error bars of r,d hsted in Tables HVl and Ivml were 
used to calculate the standard deviation of the func- 
tion. The = 17 and P^ « 96% between Eq. ^ and 
the measured R25 [not including the (rgd,i?25) — (0,0) 
point]. Since P^ is greater than the 5% significance level, 
this test does not invalidate Eq. ((921) . The R25 is a func- 
tion of Tgd and u. The r^d is a parameter of each galaxy, 
the u is a galaxy classification parameter, and the "2" 
implies the Principle of Change is applicable to the R25 
parameter. 

If e is the total energy entering a galaxy and luminosity 
(^^Q-OAMj -J ^Yie total energy emitted from a galaxy, then 
e oc 10^° "'*^J in a stable galaxy, where Mj is the total B 
band absolute magnitude of the j*'' galaxy and B band 
luminosity is a fixed proportion of the total luminosity. 
TuUy and Fisher j69l | originally suggested luminosity is 
proportional to its effective mass of the galaxy. Then the 
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FIG. 13: Plots ol the I — r profile for the sample galaxies without a rotation curve. The straight line is drawn between the 
chosen pixels. The slope on the line indicates the asymmetry. 



force Fp (ergcm^^ s^^ Mpc^^) exerted by the mass and 
p of the neighboring galaxies on the target galaxy [tJI 
along the polar axis is 

Fp^iff^X^lO-^^^.-Lz,, (93) 

where Kf^ is the proportionality constant, the sum is over 
the closest galaxies, n is the number of galaxies included 



in the calculation, Rj is the vector distance from the 
target galaxy to the j*^ close galaxy, and Zj is compo- 
nent of Rj along the target galaxy's polar axis. In the 
calculation herein, Kfz = 1 and n = 10. 

The calculation followed Hodge and Castelaz [tJI- IC 
1613, NGC 0224, and NGC 0300 were omitted from the 
calculation since the Milky Way is one of the closest 
galaxies and the M of the Milky Way is uncertain. The 
result of the calculationlssj for 28 of the sample galaxies. 



30 



I 

11 




■1 :x. 



FIG. 14: Plot ol W20 The filled diamonds use 

the measured values of W2Q and Vrmax- The line is a plot of 
Eq. 




COO 0.07 0.14 
^<zr- (l<po) 

FIG. 15: Plot of the measured R25 (kpc) versus the distance 
r^d (kpc) along the major axis of a HST image to the first dis- 
continuity or infiection point using the Dc stated in Table UTTI 
and IVllH The filled diamond, filled square, filled triangles, 
fiUed circles, open diamonds, open squares, open triangles, 
and open circle are data points along a linear R25 - r^d line 
using data from Table IDO and correspond to « = 1, it = 2, 
u = 3, u = 4, u — 5, u — 6, u = 7, and « = 8 as listed in 
Table 1X1 respectively. 



The Fp error was calculated using 0.2 mag error for Alj 
and a 5% error for distance. 

FigurelTBlis a plot of Fp versus r^d- Each hne in Fig. [TBI 
is a different "classification". The {r^d,Fp) = (0,0) point 
was included in all classifications. Table Ixll lists the cal- 
culated values for the lines of Fig. ^| 




O.UC 0,07 0,-4 
r^,.-; (Kpc J 

FIG. 16: Plot of the measured Fp (ergcm"^s~^ Mpc~^) 
versus the distance red (kpc) along the major axis of a HST 
image to the first discontinuity or infiection point using the 
Dc stated in Tables UTTl and 171111 The filled diamonds, filled 
squares, filled triangles, filled circles, open diamonds, open 
squares, and open triangle are data points along a linear Fp 
- Ted line using data from Table HXI and correspond to w = 1, 
w = 2,w = 3, w = 4, w = 5, w = 6, and to = 7 as listed in 
Table DTTI respectively. 



TABLE XI: Data of the lines in Figures [T^ 





" line' 


Corr.' 


F Test 






1 


filled diamonds 


0.99 


0.99 


13.2± 0.5 


-0.04±0.03 


2 


filled squares 


0.94 


0.89 


42 ± 3 


0.7 ±0.7 


3 


filled triangles 


0.93 


0.86 


97 ± 1 


-0.07±0.04 


4 


filled circles 


0.99 


0.99 


185 ± 6 


0.4 ±0.3 


5 


open diamonds 


0.99 


0.98 


680 ±30 


-0.5 ±0.6 


6 


open squares 


0.93 


0.91 


1570 ±90 


3 ±2 


7 


open triangle 


ND9 


ND 


4570 







a 


The integer denoting 


; the value of w 


in EQ.I95I 





'The identification of the line symbol in Fig. 1161 
'^The correlation coefficient. 

''The least squares fit slope of the lines in Fig. 1161 The units of 
Ks{ are xlO'' ergcm"^ s~^ Mpc~-^ kpc~^. 

•^The least squares fit intercept of the lines in Fig. 1161 The units 
oi Ki„ are xlO^ ergcm"^ s~^ Mpc~^. 

•^The lines are calculated using (r^di-Fp) = (0,0) as a data point. 

9 "ND" means no data because there are only two points (one data 
point). 



The lines of Fig. El are straight hnes. Therefore, the 
equation for each line is 

Fp^K,fr,d + Ki^, (94) 

where Kgf and ifjw are the slopes and intercepts of the 
lines in Fig. 1161 respectively. 

Unlike the previous r^d relationships, i^iw deviates 
slightly from zero. However, the slopes among the lines 
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in Table IxH obey the relation 



log 



10 



xlO^ ergcm"-^ Mpc kpc 



(95) 

where ifsw = 0.41 ±0.01 and Ki^ = 0.73 ±0.05 at la are 
the least squares fit of the slope and intercept of the lin- 
ear relation of Eq. (|95l) , respectively. The w is an integer 
of one to seven depending on the line from Table IXll The 
correlation coefficient and F test for Eq. (|95|l are 0.998 
and 0.996, respectively. Note, the lowest slope corre- 
sponds to = 1 in Eq. (|95|l . The value of logio(e) « 0.43 
is within 2a of Kgw 

Therefore, restating Eq. H95() as a strong Principle of 
Repetition and combining Eqs. (|94|l and (|95|) yields 



where Kp^ — 
ai la and ii'iw 



(4.5 ±0.8) X 10*^ ergcm"2 



(96) 

s"iMpc~2]jp(,-i 



0. 



The error bars of Fp listed in Tables IIXI were used 
to calculate the standard deviation of the function. 
The = 27 and w 34% between Eq. ^ and the 
measured Fp [not including the (r^d^Fp) — (0,0) point]. 
Since is greater than the 5% significance level, this 
test does not invalidate Eq. . The Fp is a function of 
Ted and w. The red is a parameter of each galaxy, the w is 
a galaxy classification parameter, and the "e" implies the 
Principle of Repetition is applicable to the Fp parameter. 

Since e = K^^ fed, 



Ted = 



(97) 



where i^ed is the proportionality constant. 

Since Zj may have any value, the continuous factors of 
the Ted relationships may be due, in part, to Fp. 

Figure 1171 is a plot of galaxy luminosity L versus r^d ■ 
The L value is calculated as 



L = 10 



-0.4 Mv 



where 



Ah 



mb -25 + 5 logio Dc 



(98) 



(99) 



where the "-25" and "5" in the equation are because the 
unit of Dc is Mpc. 

Each line in Fig. El is a different "classification". Be- 
cause six of the nine classifications had less than four 
galaxies, the {r^d,L) — (0,0) point was included in all 
classifications. Table IXIII lists the calculated values for 
the hues of Fig. [T7I 

The lines of Fig. El are straight lines. Therefore, the 
equation for each line is 



L = KsLrred + Kii^r, 



(100) 



where K^i^,- and K^^^ are the slopes and intercepts of the 
lines in Fig. 1171 respectively. 



1 




FIG. 17: Plot of the measured L (ergcm~^ s~^) versus the 
distance r^d (kpc) along the major axis of a HST image to the 
first discontinuity or inflection point using the Dc stated in 
Tables Hm and rTim The filled diamond, filled squares, filled 
triangles, filled circles, open diamonds, open squares, open 
triangles, open circles, and "X" are data points along a linear 
L - r^d line using data from Table llXl and correspond to x = 1, 
X = 2, X — 3, X = 4, X = 5, X = 6, X — 7, X — S, and a; = 9 as 
listed in Table IXTTl respectively. 



TABLE XII: Data of the lines in Figures [T7I 



x" 


line' 


Corr.'^ 


F Test 






1 


filled diamond 


ND9 


ND 


26 


0.0 


2 


filled squares 


0.99-1- 


0.99+ 


48± 


0.0+0.0 


3 


filled triangles 


0.994- 


0.99+ 


68± 1 


-0.1+0.2 


4 


filled circles 


0.99-1- 


0.99+ 


133+ 5 


-0.1+0.2 


5 


open diamonds 


0.97 


0.92 


285+ 20 


0.1+0.7 


6 


open squares 


0.994- 


0.99+ 


474+ 18 


0.2+0.4 


7 


open triangles 


0.97 


0.94 


851+100 


2 +2 


8 


open circles 


0.99-1- 


0.99+ 


2100+100 


-0.2+0.3 


9 


X 


0.99+ 


0.99+ 


4000+200 


0.3+0.6 



"The integer denoting the place of the line in the order of increas- 
ing ifsLr- 

''The identification of the line symbol in Figures ^lYI 



■ r^tj lines. 

■ r^d lines. 



"^The correlation coefficient. 

''The least squares fit slope (ergcm~^ s"'^ kpc) of the L - 
''The least squares fit intercept (ergcm~^ of the L - 
•'The lines are calculated using (r^jjii) = (0,0) as a data point. 
9 "ND" means no data because there are only two points (one data 
point). 



Since i^iLr ~ 0, the slopes among the lines in TablePnTI 
obey the relation 



log 



KsLr 



10 I _9 

' erg cm ! 



"1 kpc 



Kf,^x + Kh 



(101) 



where if^x = 0.274+0.007 and Ki^ = 1.08±0.04 at la are 
the least squares fit of the slope and intercept of the linear 
relation of Eq. (|101|l . respectively. The x is an integer of 
one to nine depending on the line from Table Ixlll The 
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correlation coefficient and F test for Eq. (|101(l are 0.99 
and 0.99, respectively. Note, the lowest slope corresponds 
to x = 1 in Eq. H101|l . The L/r^d is not independent of 
the distance used to calculate the L and r^^. The value 
of logio(2)Ri 0.301 is within 4cr of Kgyi- There are only 
three galaxies on the x ^ 1 and x = 2 lines. If the x = 1 
and X = 2 lines are ignored, the K^x — 0.29 ± 0.03 that 
is within la of logio(2). 

Therefore, restating Eq. IjlOlfl as a strong Principle of 
Change and combining Eqs. (|100|l and (|101|l yields 



value of Ki = Y^/{2Xi) with an Xi > 0.4. (5) Calculate 
the ti value for each set of (Xi,Yi) values as 



L = i^Lr 2" r,d, 



(102) 



where Ki,i- — (8.9 ± 2.0) erg cm ^ s ^ kpc at la and 
J^'ir25 = 0. Note the lack of data points in the upper 
right of Fig. [13 

The error bars of r,d hsted in Tables Hvl and IVlTIl were 
used to calculate the standard deviation of the func- 
tion. The = 30 and ss 44% between Eq. ((TU^ 
and the measured L [not including the [r^d, L) = (0,0) 
point]. NGC 4639 had a = 18. Since is greater 
than the 5% significance level, this test does not invali- 
date Eq. (|l()2|l . The L is a function of r^d and x. The 
Ted is a parameter of each galaxy, the x is a galaxy clas- 
sification parameter, and the "2" implies the Principle of 
Change is applicable to the L parameter. 

The Fp is calculated from the magnitude, which is de- 
rived from e, and the vector position of other Sources. 
Since Fp and M are independent variables, r^d of a galaxy 
depends on the Mj, hence e, of all galaxies. Also, other 
galaxy parameters are related to r^d. This suggests the 
e and the position in Euclidian space of all the Sources 
(center of galaxies) completely determines the parame- 
ters of a target galaxy. 



V. THE MANY LINES IN THE PARAMETER 
RELATIONS WITH r,d ARE NOT RANDOM. 

The many lines in each plot in Figs. IBl [TUl [TTl [HI [T5l 

1161 and 1 171 mav suggest the relationships portrayed may 
be just random data points. The null hypothesis to be 
tested is "the data points are indeed random points" . 
This null hypothesis may be tested by following the pro- 
cedure used in discovering the Eqs. ([S3|), JTTI, 113, l|7^ - 
([^ . and ([TU^ as follows: (1) Generate 15 sets of 

two random numbers between zero and one. Call one of 
the i^^ set the independent variable {Xi) and the other 
the dependant variable (Yi), where i varies from one to 
15. (2) The equation to be tested is 



Yr 



calc 



(103) 



where Ki is the proportionality constant and ti is the 
integer classification for the i^^ set. (3) A trial consists 
of generating the random numbers for the 15 sets and 
subjecting the sets to the slope, correlation coefficient, 
and calculations and tests the galaxy samples passed. 
(4) Calculate the Ki value for each trial as the minimum 



t, = ROUND 



log2 



Y 



kiX, 



(104) 



where "ROUND" means to round the value in the braces 
to the nearest integer. (6) If any one of the ti > 7, the 
trial failed. (7) For calculating the slopes and correlation 
coefficients, if any subset lacks a sufficient number of val- 
ues, that subset will be ignored. (8) If the correlation 
coefficient of any ti subset of {Xi,Yi) values, including 
the {Xq^Yq) = (0,0) point, is less than 0.95, the trial 
failed. (9) If the linear slope of the relation between the 
ti and slope, including the [Xq, Yq) = (0,0) point, of the ti 
subsets 7^ (0.30 ±0.03), the trial failed. (10) If the corre- 
lation coefficient of the line in (9) above < 0.95, the trial 
failed. (11) Calculate the Fcaici according to Eq. (|103|l . 
If the P^ < 0.05 between the Yi and Fcaici for the 15 sets, 
the trail failed. (12) maintain a count of the number N^^ 
of sets of {Xi,Yi) with Xf + Y^^ < 1 and the number 
TVxitY of sets of {X„Yi) with X, < Y,. (13) Redo the 
trial with another 15 sets of random numbers for 30,000 
trials. 

The results are: (1) The iV^/(15x30, 000) = 0.78528 « 
7r/4 and 7VxitY/(15 x 30,000) = 0.49999. Therefore, the 
random number generator performed satisfactorily. (2) 
22 of the 30,000 trials passed (0.07%), the remainder of 
the trials failed. Therefore, the null hypothesis is rejected 
with a 99+% confidence. 



VI. DISCUSSION 

Space has a physical presence, a density, and an ability 
to exert a force on the a of hods (matter) without be- 
ing matter. The Space in this Paper looks like the ether 
with added properties. The Michelson-Morley experi- 
ments "t^ applied to the velocity of photons. To achieve 
the results of Michelson-Morley's experiments, photons 
must have zero surface area presented to their direction 
of travel. Hence, photons must be hods arranged in a col- 
umn. However, recent re-examination of the Michelson- 
Morley data [73 and more accurate measurements [t^ 
suggests a preferred reference frame may exist. The CUM 
suggests these findings indicate ripples on the hods of 
photons provide a dragging force on the photon. There- 
fore, in addition to testing the ether model, Lorentzian 
relativity, and Einstein's special relativity, the Michelson- 
Morley type experiment may be testing the p, Vp, and 
d(Vp)/di in our solar system. The ratio of the diameter 
of the earth's orbit to the distance from the center of the 
Galaxy suggests the Michelson-Morley experiment may 
require greater precision and appropriate timing to test 
the p, Vp and d(Vp)/di. 

Ferrarese ,39] discovered a linear relationship between 
circular velocity Vc beyond R25 and bulg e velocity disper- 
sion (Tc (see her Fig. 1). Baes et al. [23 expanded on the 
data. NGC 0598 was a clear outlier. Galaxies with a^ 
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less than 70 kms"""^ {vc less than 150 kms~^) also deviate 
from the linear relation. Also, galaxies with significant 
non-circular motion in the HI gas such as NGC 3031 were 
omitted in Ferrarese The Vc for NGC 0598 used in 
Ferrarese 3.9j was 135 kms^^, which is the highest data 
point of a rising RC ■ Most of the other curves Fer- 
rarese used were flat. For flat RCs Vc « Vrmax- If 
the Vc = Vrmax = 85 kms^^ is used by Ferrarese [s^ 
for NGC 0598, NGC 0598 will fit on her plotted line. 
Thus, the range of v^ — relation may extend to 27 
kms-i < CTc < 350 kms"^ NGC 3031 has a declin- 
ing RC, i?25 ~ 11 arcmin (from the LEDA database), 
VMI25) « 190 kms-i and CTc = 174 ± 17 kms'^ 
1^. NGC 3031 is not on the plotted line of Ferrarese 
[29'|. However, if Vc — Vrmax = 243 ± 14 kms~^ is used, 
NGC 3031 does fit on the plotted line. 

Multiple slopes of the Ca versus Ep, (erg cm"'^ s~^) 
relation from Hodge and Castelaz 'TlJ (see Hodge and 
Castelaz 71] for the definition and calculation of Ca and 
Es,) also obey Eq. Separate the CETl zone into two 
zones. The first zone (CETla) has NGC 3621 and NGC 
5253 which have low E^, and a slope of — 2.28 x 10^^ erg^^ 
cm2 s\ The second zone (GET lb) has NGC 3031 and 
NCG 2090 which have higher E^ and a slope of 8.87 x 
10"'' erg~^ cm^ s^. The slope values were arranged in 
ascending order with CETlb, CET2, GETS, and CET4 
corresponding to n=l, n=2, n=3, and n=4, respectively. 
For the Ca versus i?a relationship. Unit is erg~^ cm-^ s^, 
Ksn = 0.585 ±0.005 and Kin = -0.54 ±0.02 at 1 a. The 
correlation coefficient and F test for Eq. (|52|l are 0.99 and 
0.99, respectively. 

The linear potential model (LPM) explored by 
Mannheim and Kmetko [b^ and references therein] is 
derived from relativistic considerations. The LPM pro- 
poses a gravitational potential V{r) of the form V{r) = 
—P/r+jr/2 where /3 and 7 are constants. The potentials 
in the CUM are proportional to r^^. RCs with u^j.^. oc rs„ 
or Vfj.j. oc r|.j. are derived by geometric and particle type 
variations after the force of the Mkmax is balanced by the 
Space force. However, both models conclude the intrinsic 
Vsrr and Vs are a linear function of r. Hence, both pre- 
dict rising rotation curves. The LPM derives flat curves 
by mass distribution but ultimately has rising curves. If 
LPM is to derive declining curves, it is by the influence of 
other galaxies which have a linear r potential. The LPM 
does not predict the relations with red or oc r^.^.. 

The DM model (DMM) places large amounts of un- 
known, non-luminous matter in a halo around a galaxy 
to explain rising rotation curves. Whereas, the DMM 
requires the DM term to become larger for rising rota- 
tion curves, CUM requires the r • term to become 
smaller. The disagreement between current CDM the- 
ory and observation is most evident in galaxies with low 
surface brightness (LSB) such as NGC 0300 and rising 
rotation curves such as NGC 0300, NGC 0598, and NGC 
3319. These galaxies fit in the CUM. In addition, the 
DMM would consider Urmax as merely another point on 
the rotation curve. The DMM suggests mass is infalling 



from outside the galaxy. The CUM suggests matter is 
outflowing as gas and inflows as denser particles. The 
DMM is not really a predictive model but rather a pa- 
rameterization of the difference between observation and 
Newtonian expectation. Because of the number of pa- 
rameters in the DMM, it is very difficult to falsify. How- 
ever, the large number of parameters makes the DMM 
cumbersome and, therefore, replaceable. 

Padmanabhan and Choudhury fldi and references 
therein] explored the possibility of a scalar field acting 
as dark matter. The CUM suggests the massless (mass 
density = 0) scalar p field at small scales causes an intrin- 
sic rising RC. The DMM would require a large amount 
of DM. The p fleld acts like large amounts of DM at 
galaxy diameter scales. Therefore, the CUM suggests 
the amount of DM in the universe is zero. Also, this 
Paper suggest the p field causes Newtonian kinematic 
measurements to considerably underestimate mass in a 
galaxy. For example, Eq. with Mcs = M, suggests 
the larger Mtckmax is the more accurate measure of the 
amount of mass in the central region of galaxies. 

The CUM's {rUgMg — mse)/r term has the form of a 
subtractive change in acceleration. The CUM suggests 
this causes the species of matter to change with radius 
which causes the appearance of a gravitational acceler- 
ation change with radius. In the MOND model 0, the 
Newtonian gravitational acceleration is measured in our 
solar system, which is in the low r end of the SR. The 
MOND data comparison is usually in the far RR and 
SR regions. Hence, the existence of a differing e/r or 
an apparent modified gravitational acceleration is ex- 
pected by the CUM model. NGC2841 is a candidate 
to falsify MOND. A distance greater than 19 Mpc, com- 
pared to the measured of 14.1 Mpc |^, is needed for 
MOND 0. However, NGC2841 data fits the CUM. Fur- 
ther, given a relatively short radius range of the MOND 
application, the e/r change may be approximated as a 
constant. Also, since we are near the Wrmax in the Galaxy 
and since the CUM holds the relationship of v'^^^^ — e 
as a universal relationship, the MOND proposition that 
the modified gravitational acceleration is a universal con- 
stant follows. If in the SR there is another acceleration 
from r • Oo to balance the e decline, then the MOND's 
fiat rotation curve can result. MOND also predicts that 
if the rotation curve is flat, it is axis symmetric. Flat 
and declining rotation curves are asymmetric and 
references therein] . Rising rotation curves are more sym- 
metric. Therefore, the CUM reduces to MOND where 
Ar • tto ^ and V^{r • Fo) = in the SR. 

The model that places a SBH at the center of galaxies 
says nothing about the RC outside the KR. Also, this 
model of the center of galaxies seems to be invalidated 
by observations of gas near the center of the Galaxy 129, 
page 595]_y0| i by the apparent inactivity of the SBH [S^, 
133. 133I l34| , by the observed relations with other galaxy 
parameters [S 111 111 IH H El El EHii El , andby 
the apparent need for a feedback control mechanism ^5] . 

The DMM adds several parameters per galaxy. Both 
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the LPM and MOND add one parameter to Newto- 
nian mechanics and another mildly varying parameter 
to achieve their fits and have been applied to only close 
galaxies. The CUM adds one parameter (e) per galaxy 
and one term per nis species to Newtonian dynamics, 
but needs other galaxies and their distances to further ex- 
plain a galaxy's dynamics. The e of other galaxies modify 
the equations to allow flat and declining RCs and to al- 
low asymmetries in RCs. Neither the LPM nor MOND 
address the relationships outside the RR and SR. The 
DMM, LPM, MOND, and SBH galaxy models either are 
invalidated by or fail to explain the data examined herein. 

I consider it very attractive that the CUM, at this 
early stage of model development, has brought together 
heretofore apparently unrelated observations and has dis- 
covered new galaxy parameter relations. 



VII. CONCLUSION 

A changing universe model (CUM) was developed in 
accordance with a new set of Fundamental Principles. 
The CUM posits the stuff of our universe is continually 
erupting into our universe from sources at the center of 
galaxies. The stuff of our universe is matter and a mat- 
terless scalar whose gradient exerts a force on the cross 
section of matter. An energy continuity equation (Sec- 
tion was derived. Equations describing the forces 
of gravity fSection lllB 3|) and Space fSection lll B l|l were 
derived. A description of the structure of particles (Sec- 
tion was presented to provide an explanation of the 
phenomena observed at the Galaxy center. An expla- 
nation of Galaxy structure was presented. The CUM 
qualitatively suggests: (1) That our universe is flat. (2) 
The characteristic of the right angle is orthogonality. (3) 
That total kinetic energy equals total potential energy in 
our universe. (4) That Vm/C, a, t, and l are constants of 
the universe. (5) That Space is causally connected and 
coherent. (6) Why there are three, and only three. Space 
dimensions. (7) Why there are three, and only three, 
families of particles. 

To test the proposition, a single equation was derived 
using the CUM that calculates the rotation velocity v of 
particles along the major axis of spiral galaxies. The 
equation adds one parameter for each mature galaxy 
and one parameter for each particle (atom nucleus and 
smaller) species to Newtonian dynamics. The equation 
is used to correlate the measured mass M, of the theo- 
rized, central supermassive black hole (SBH) with other 
galaxy parameters and to describe the rotation velocity 
curve (RC) of spiral galaxies. The application of the 
CUM suggests galaxies are divided into regions. The 
CUM suggests the intrinsic rotation curve in the disk re- 
gion is rising and suggests flat and declining disk rotation 
curves are cause by neighboring galaxies. Each region of a 
galaxy has its own physics. The single equation, through 
the relative value of terms, reduces to a description of the 
physics of the RC of each region. Four new parameters 



were defined: (1) the radius r^d to a discontinuity in the 
surface brightness versus radius curve, (2) the rotation 
velocity "yrmax at the maximum extent of the rapidly ris- 
ing region of the RC, (3) the radius rmiax to the Wrmax, 
and (4) the Source strength e of each galaxy. 

In a sample of 20 galaxies with the radius of influence 
of the central black hole resolved, this Paper predicted 
a linear Af, — r^d relation and found A/, — A'sc?'ede'^, 
where AT^o = (3.7tj]:J^) x lO^Mgkpc^i, e « 2.718 (the 
natural log base), and n is an integer from one to six 
depending on the galaxy. Since the sample included E, 
SO, and S galaxies, the derived equation blends smoothly 
across the galaxy morphologies of the sample. 

In a sample of 15 spiral galaxies with Hi RCs and 
Cepheid based distances, this Paper predicted a lin- 
ear Trmax Wmiax " ^ed relation and found rj-max IJrmax = 

Ksorr.d e^, where AT^er = (4.8 ± 0.9) x 10^ km^s'^ and 
m is an integer from one to six depending on the galaxy. 
Also, this Paper predicted a linear rrmax — ''ed relation 
and found r-rmax = AT^ordt ''ed 2P, where ATscrdt = 36 ± 8 
and p is an integer from one to seven depending on the 
galaxy. Also, this Paper predicted a linear Wrmax ~ ''ed re- 
lation and found l^rmax = A'sorvds »'ed 2^ whcrC A'sorvds = 

(1.6 ± 0.3) X 10^ km^ s^^ kpc~^ and s is an integer from 
one to seven depending on the galaxy. Further, for each of 
the four linear relationships, linear relationships between 
the respective parameter and r^d for each of the integer 
classifications were also found. The CUM predicts the 
square of the rotation velocity is linearly related to r in a 
portion of the RR and is linearly related to in another 
portion of the RR. 

For a sample of 31 spiral galaxies with a wide range of 
characteristics, the Hi (21 cm) line width W20 (kms~^) at 
20% of the peak value corrected for inclination has been 
found to be tightly related to Urmax by lU2o/kms^^ = 
ATi (^>rmax/kms~l)^^ where ATi = 7.3, K2 ^ 0.78, and 
the uncertainty in W20 is ±12% with « 86%. Since 
TU20 is related to the total absolute magnitude Mg (mag.) 
in band f3 by the TuUy-Fisher relationship ^, 70], 
is also related to r^d- The de Vaucouleurs radius R25 
(kpc) is proportional to r^d by R25 = ATpir 2" r^^ with 
« 96%, where ATpr = (32 ± 7) at la and u is an in- 
teger from one to eight depending on the galaxy. For 28 
of the 31 sample galaxies, the force Fp of the e, which 
is proportional to the luminosity, of neighboring galaxies 
directed along the polar axis of the target galaxy is re- 
lated to Ted by Fp = ATpre™ r^d with P^ « 34%, where 
A'pr = (4.5 ±0.8) X 10^ ergcm-^s-iMpc-^kpc-i at la 
and w is an integer from one to seven depending on the 
galaxy. This suggests the continuous factors of the r^d 
relationships are due, in part, to Fp. For the 31 sample 
galaxies, the luminosity L in the visible band (from the 
NED database) is proportional to r^d by L = ATLr 2^ r^d 
with P^ « 44%, where ATLr = (8.9±2.0) erg cm'^ s"^ kpc 
at la and x is an integer from one to nine depending on 
the galaxy. 

The Ted is related to galaxy dynamics parameters, lu- 
minosity parameters, and an isophotal radius parameter 
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by the Principle of Change. Also, r^d is related to Fp and 
effective mass parameters by the Principle of Repetition. 

The CUM is not invalidated by the data presented 
herein. 
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APPENDIX A: PRINCIPLES 

The Principal of Fundamental Principles is that a Fun- 
damental Principle, and the models developed from it, is 
a meaningful and useful principle that applies to all levels 
of physical systems. To be meaningful is to be able to be 
used to describe and predict the outcome of experiments. 
To be useful is to be able to be used to cause desired out- 
comes. An outcome of an experiment includes placing 
bounds on what (parameter or event), where (position 
coordinates), and when (time coordinate). The level of a 
physical system refers to the size of the domain of appli- 
cability over which a set of physical theories applies such 
as galaxies versus atoms and Newtonian versus special 
relativity. 

Corollary I is if a candidate to be a Fundamental Prin- 
ciple is found to not apply in some system then it is not 
a Fundamental Principle. Corollary II is if a principle is 
found in all of physical systems, then it is likely to apply 
to larger and smaller levels and to new concepts. 

The Principle of Superposition, the Correspondence 
Principle, and Principle of Minimum Potential Energy 
are such Fundamental Principles. The Correspondence 
Principle is an interpolation of the Fundamental Princi- 
ples. Corollary II is an extrapolation of the Fundamental 
Principle. 



A "scientific model" (theory) is derived from the tran- 
scendent idea of the Fundamental Principle and is appli- 
cable to a defined domain. Since life and social systems 
are physical systems, by Corollary II and the Cosmo- 
logical Principle, the transcendent idea of the Funda- 
mental Principles must also apply to life and to social 
systems. The more fundamental scientific models have 
larger domains. A proposal becomes a scientific model 
when a deduction from it can be found to be factual (is 
observed). This concept does not require a candidate 
model to include "falsifiable predictions" and does not 
invalidate a scientific model because of the existence of 
falsifying observations. For instance, Newtonian dynam- 
ics is a valid scientific model. Observations in the do- 
main including relative velocities approaching the speed 
of light falsify Newtonian dynamics. However, this only 
limits the Newtonian dynamics domain. Religious ideol- 
ogy, models based on belief, and philosophy models may 
be scientific models provided they are useful and mean- 
ingful with restricted what, where, and when bounds. To 
survive, a scientific model must compete for attention. 
The concept of a scientific model survives because the 
human mind is limited in the ability to maintain a cata- 
log of the nearly infinite number of possible observation. 
Scientific models with empty or more limited domains 
have little usefulness and meaningfulness. 

The Copernican Principle, which states we are not 
privileged observers of the universe, is accepted. In quan- 
tum mechanics our presence may change the outcome of 
experiments. However, this is true for any observer. If, in 
some set of observations, we appear privileged, the priv- 
ilege must be incorporated in the model. For example, 
we are in a galaxy disk and are close to a sun. We are in 
a highly unique and privileged area. Just because we are 
carbon based does not imply all intelligent life is carbon 
based. 

The Anthropic Principle is accepted to the extent 
that what is observed must have been created and have 
evolved to the present. What is observed must be able 
to be observed. Note this statement of the Anthropic 
Principle omits the requirement that it depend on an 
assumption of "life" and "intelligence" since life and in- 
telligence are inadequately defined. The existence of life, 
social systems, and intelligence are observations of our 
universe and, therefore, must be able to exist. An unob- 
served parameter may or may not be able to be observed. 
Therefore, the negative model candidates are not useful. 

Some argue |77l and references therein] the Anthropic 
Principle cannot be part of science since it cannot yield 
falsifiable predictions. In this statement of Fundamental 
Principles, what is observed can be used to define the 
domain of scientific models. 

The Cosmological Principle is accepted to the extent 
that the physics must be the same at all positions in the 
universe. Therefore, the theories must explain the con- 
ditions at all positions in the universe. For example, the 
physics that states all objects fall to earth was found to 
be limited when Galileo noted moons rotating around 
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Jupiter. Newton's restatement of physics was more cos- 
mological. However, the Cosmological Principle is not 
extended to imply the universe is isotropic and homo- 
geneous. Near a star is a unique position. The physics 
must explain other positions. The physical theories must 
explain any isotropics and anisotropics. 

The Principle of Change states that each system is 
built of previous systems by a minimum step change. 
When the size of a structure of an object becomes lim- 
ited, a new structure comprising a combination of exist- 
ing elements can occur. Change is a building on existing 
structures. Destruction of objects into its parts docs not 
foster the Principle of Change. Since the minimum step 
Change is a split from one to two, a corollary (strong 
Principle of Change) is that whatever exists will split 
into two parts. The two parts may be a continuation 
of what was plus a new system or may be what was in- 
corporated into the two new systems. Particles became 
a hydrogen atom followed by evolution of other atoms. 
Atoms became molecules. A model that requires a large 
step where there are possible intervening steps is not ob- 
served and is forbidden. 

A corollary of the Principle of Change is that all things 
have a beginning. 

The Principle of Repetition states that there are two 
ways to repeat a Change: (1) If conditions allow an ob- 
servable change, then the change will occur again un- 
der similar conditions. (2) The repeated Changes have 
a common cause (reproduction). A corollary is that if 
two systems have the same observable results, then sim- 
ilar conditions exist or the systems were reproduced. A 
strong statement of the Principle of Repetition is that 
the amount of increase of a parameter by the Principle 
of Repetition depends on the size of the parameter. De- 
struction of objects to have "room" for "the new" is a 
Repetition because the only objects that can be built 
from the pieces are a Repetition of objects. 

The Principle of Negative Feedback states that any 
system with a relatively narrow parameter relationships 
must evolve from a broader system and must have a neg- 
ative feedback loop to maintain the narrow parameters 
and achieve balance between the Change and Repetition 
processes. Otherwise, the system is unstable and tran- 
sitory. The effects of the unstable system will cease to 
exist without consequential fallout or permanent change. 
We observe objects that have limited size. So, there is a 
limiting force or negative feedback condition controlling 
the size of ciach object. So too must black holes have a 
size limitation and a negative feedback condition. 

The Principle of Minimum Potential Energy can be 
stated as a Principle of Minimum Action, which states 
the path of least energy expenditure will be followed dur- 
ing the change from one state to another. Combining this 
with the Principle of Negative Feedback yields the con- 
clusion that the feedback energy equations are linear, of 
the first time derivative, and of the first time integral. 
In an equation none, one, or two of these elements may 
be zero. Since a higher power feedback equation would 



result in the same controlled output and consume more 
energy, the minimum equations are sufficient. Therefore, 
the assumption of "linear" in formulating descriptive, 
negative feedback energy equations is sufficient. 

The Principle of Geometric Rules states that the ob- 
served geometric relationships apply in all levels of sys- 
tems. Hence, the conservation of energy/mass must be 
related to geometric rules we observe in our universe. 
Hence, tt = circumference / radius in two dimensions 
must be the same number in three dimensions, tt is not 
only a universal constant in our universe, it is a constant 
in all the universes. The division by two is another uni- 
versal concept. The division by two for each dimension 
into equal angles yields the right angle. 



APPENDIX B: VISUAL BASIC 5.0 PROGRAM 
FOR FINDING ra. 

This program starts with the appropriate information 

extracted from the HST FITS file and the operator has 
input the required additional data. 

Option Explicit 

Const Pi As Single = 3.141592654 

'7. integer 

'& long integer 

' ! single precision 

'# double precision 

'$ string 

Dim it, j7., x%, y%, k&, 1&, m&, p7. 
Dim d!(802, 802), dmax ! , Xmax!, Ymax! 
Dim XL!, YL!, RAp#, DECp#, Dist#, PA! 
Dim Ap#, bO#, cp#, Xc#, DECO# 
Dim XLvfe, YLvfe, XLvLast7., YLvLast7., RAO# 
Dim crVall#, crVal2#, crPixl#, crPix2# 
Dim NAXIS17., NAXIS27., ORIENTAT!, CD1_1# 
Dim dmaxi!, XmaxI ! , YmaxI!, CD1_2#, CD2_1# 
Dim VegaZP#, ExpTime#, Area*, CD2_2# 
Dim Instrument!, Filter$, IData! (2270) 
Dim AData#(2270) , Sb!, TSlope!, Rip! 
Dim SbO!, CountL7., CountH7., Rln! , Real! 
Dim Rea2!, Sxy#, Sx#, Sx2#, Sy# 
'crVall, crVal2, crPixl, crPix2, ExpTime 
'Instrument, Filter, 

'NAXISl, NAXIS2, ORIENTAT, CD1_1, CD1_2 
'CD2_1, CD2_2 

'are as defined in the HST FITSFILE header. 
'VegaZP is the zero point of the surface 
'brightness calculation, input by operator 
'd(X,Y) is the data number (DN) 
'of the pixel 

'at the X Eind Y image coordinate. 
'Since the galaxy may be asymmetric, 
'd(X,Y) is not averaged or 
're-calculated statistically. 

Private Sub cmdDoOne_Click() 



PA = Clnt(Textl.Text) 'position 
' angle of galaxy, operator input 

Area = Abs(CDl_l * CD1_2) * 3600 _ 
* 3600 

'find center of gal. by highest pixel 
' value close to input coordinates 
Xmax = Val(Text2.Text) 
Ymax = Val(Text3.Text) 
dmax = dCXmeix, Ymax) 
DoAgain2 : 

For X = Xmax - 20 To Xmax + 20 
For y = Ymax - 20 To Ymax + 20 
If X > And y > And x < _ 
MAXISl And y < NAXIS2 Then 

If d(x, y) > dmax Then 
dmax = d(x, y) 
Xmeix = x 
Ymeix = y 
GoTo DoAgain2 
End If 
End If 

Next 

Next 

Text2.Text = CStr(Xmax) 

Text3.Text = CStr(Ymax) 

PA = PA - ORIENTAT 
' PA is now angle from X axis 
'counterclockwise on image and Y axis 
'is up from bottom 
DoAgain: 

If PA >= 180 Then 'not matter 
' for tcin function 

PA = PA - 180 
GoTo DoAgain 
End If 

If PA < Then 

PA = PA + 180 

GoTo DoAgain 
End If 

If PA > 89.3 And PA < 90 Then 
'TAN (90) is problem 
PA = 89.3 
End If 

If PA < 90.7 And PA > 90 Then 
PA = 90.7 

End If 

'record general information in Datal 
' which is an Excel file 

Datal .Refresh 

Datal . Recordset . MoveFirst 

Datal . Recordset . Edit 
Datal . Recordset . Fields ( "NAME" ) . Value. 
= Left(Filel.f ileneime, 
Len(Filel .filename) - 4) 
Datal . Recordset . Update 
Datal . Recordset . MoveNext 
Datal . Recordset . Edit 



Datal .Recordset .Fields ("NAME") .Value 
= Label2. Caption + " deg." 

Datal .Recordset .Update 

Datal . Recordset . MoveNext 

Datal . Recordset . Edit 

Datal .Recordset .Fields ("NAME") .Value. 
= Label3. Caption + " deg" 

Datal . Recordset . Update 

Datal . Recordset . MoveNext 

Datal .Recordset . Edit 

Datal . Recordset . Fields ( "NAME" ) . Value. 
= "IMAGE PA: " + Format(PA, "##0.#") 

Datal .Recordset .Update 

Datal .Recordset .MoveNext 

Datal . Recordset . Edit 

Datal . Recordset . Fields ( "NAME" ) . Value 
= "CNT: " + Format (d(CInt(Text2. Text) , 
CInt (Texts. Text) ) , "######") 

Datal .Recordset .Update 

Datal . Recordset . MoveNext 

Datal . Recordset . Edit 

Datal . Recordset . Fields ( "NAME" ) _ 
.Value 

= "CNT X: " + Format (Val(Text2. Text) _ 
, "###") 

Datal .Recordset .Update 

Datal . Recordset . MoveNext 

Datal . Recordset . Edit 

Datal .Recordset .Fields ("NAME") _ 
.Value _ 

= "CNT Y: " + Format (Val (Texts. Text )_ 
, "###") 

Datal .Recordset .Update 

Datal . Recordset . MoveNext 

Datal .Recordset . Edit 

Datal . Recordset . Fields ( "NAME" ) . Value 
= "INSTRUMENT: " + Instrument 

Datal . Recordset . Update 

Datal . Recordset . MoveNext 

Datal .Recordset .Edit 

Datal . Recordset . Fields ( "NAME" ) . Value. 
= "FILTER: " + Filter 

Datal .Recordset .Update 

Datal .Recordset .MoveNext 

Datal . Recordset . Edit 

Datal . Recordset . Fields ( "NAME" ) . Value 
= "VegaZP:" + CStr(VegaZP) 

Datal .Recordset .Update 

Datal . Recordset . MoveNext 

Datal . Recordset . Edit 

Datal . Recordset . Fields ( "NAME" ) . Value 
= "EXPOSURE TIME:" + CStr(ExpTime)_ 

Datal . Recordset .Update 
' initialize 

Datal .Recordset .MoveFirst 

For i = 1 To 2270 
IData(i) = 
AData(i) = 



Datal . Recordset . Edit 

Dat al . Recordset . Fields ( " INDEX " ) _ 
.Value = 

Datal . Recordset . Fields ( "X"_ 
) .Value = 

Datal . Recordset . Fields ( " Y" ) _ 
.Value = 

Datal . Recordset . Fields ( "DEGREES" ) _ 
.Value = 

Datal . Recordset . Fields ( "VALUE" ) _ 
.Value = 

Datal . Recordset . Update 

Datal .Recordset .MoveNext 

Next 

Datal . Recordset . MoveFirst 
'coordinates of gal. center 

DECO = crVall + ( (Val (Text2 . Text) 

- crPixl) * CD1_1) + ((Val (Texts. Text )_ 

- crPix2) * CD1_2) 

RAO = crVal2 + ( (Val(Text2 . Text) _ 

- crPixl) * CD2_1) + ( (Val (Text3. Text) _ 

- crPix2) * CD2_2) 

'a change of location of center on 
'the image 

If Val (Text 2. Text) > Then 
crPixl = Val (Text2. Text) 

End If 

If Val (Texts. Text) > Then 
crPix2 = Val (Texts. Text) 
End If 
p = 

If PA < 45 Or PA > 135 Then 
For XL = 1 To NAXISl - 1 _ 
Step 0.25 

XLv = Int(XL) 

YLv = crPix2 + ((XL - crPixl) _ 

* (Tan(PA * Pi / 180))) 

If YLv > 1 And YLv < NAX1S2 _ 

- 1 Then 

Record 
End If 

Next 

Else 

For YL = 1 To NAXIS2 - 1 _ 
Step 0.25 

YLv = Int(YL) 

XLv = crPixl + ((YL - crPix2)_ 

* (Tan((90 - PA) * Pi / 180))) 

If XLv > 1 And XLv < NAX1S1_ 

- 1 Then 

Record 
End If 

Next 

End If 
'record last record 

Datal . Recordset . Edit 

Datal . Recordset . Fields ( "X" ) . Value_ 
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= CSng (NAXISl) 

Datal . Recordset . Fields ( " Y" ) . Value_ 
= CSng(NAXIS2) 

Datal . Recordset . Fields ( "VALUE" ) _ 
.Value = 

Datal . Recordset . Fields ( "DEGREES" ) _ 
.Value = 

Datal .Recordset .Update 
Datal . Recordset . Close 

Re a 

Label4. Caption = "DONE: " 
+ Filel . filename 
End Sub 

Sub RecordO 

If XLv = XLvLast And YLv 
= YLvLast Then 'prevent 
'duplication of pixels 
GoTo NotDo 

End If 

XLvLast = XLv 
YLvLast = YLv 
p = p + 1 

IData(p) = d(XLv, YLv) 
'calculate the Right Assension and 
' DEClination of the pixel 

RAp = RAO + ((crPixl - XLv) * CD1_1_ 
) + ((crPix2 - YLv) * CD1_2) 

DECp = DECO + ((crPixl - XLv) 

* CD2_1) + ((crPix2 - YLv) * CD2_2) 

Ap = RAp - RAO 'left in degrees 
' for troubleshoot ease 

bO = 90 - DECO 

cp = 90 - DECp 
'calculate the cos (angle from the 
'center to the pixel) 

Dist = (Sin(b0 * Pi / 180) * 
Sin(cp * Pi / 180) * Cos(Ap * Pi / 180) )_ 
+ (Cos(bO * Pi / 180) _ 

* Cos(cp * Pi / 180)) 

'calculate arccos of Dist 

If Abs(Dist) > 0.99999999999995_ 
Then 

Dist = 

Else 

AData(p) = (Atn(-Dist / 
Sqr(-Dist * Dist + 1)) + 2 * Atn(l))_ 
* 180 / Pi 
End If 

'dist >0 unless changed 

If XLv < crPixl And YLv <=_ 
crPix2 Then_ 

AData(p) = -AData(p) 
End If 

If XLv >= crPixl And YLv <_ 
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crPix2 Tlien_ 

AData(p) = -AData(p) 
End If 

If PA > 90 Then 

AData(p) = -AData(p) 

End If 
'record DN pixel values 

Datal . Recordset . Edit 

Datal . Recordset . Fields ( " INDEX" ) _ 
.Value = CSng(p) 

Datal . Recordset . Fields ( "X" ) _ 
.Value = CSng(XLv) 

Datal . Recordset . Fields ( " Y" ) _ 
.Value = CSng(YLv) 

Datal . Recordset . Fields ( "VALUE" ) _ 
.Value = IData(p) 

Datal . Recordset . Fields ( "DEGREES " ) _ 
.Value = AData(p) 

Datal . Recordset . Update 

Datal . Recordset . MoveMext 
NotDo: 
End Sub 

Function Slope2(A As Long) As Double 
If A = Then 
Sxy = 
Sx = 
Sy = 
Sx2 = 
1 = 

Else 

Sxy = Sxy + (IData(A) 
* AData(A)) 

Sx = Sx + AData(A) 
Sy = Sy + IData(A) 
Sx2 = Sx2 + (ADataCm + k) 
* ADataCm + k)) 
1 = 1 + 1 

If Abs((l * Sx2) - (Sx * Sx))_ 
< 0.000000001 Then 

Slope2 = 

Else 

Slope2 = ((1 * Sxy) _ 
- (Sx * Sy)) / ((1 * Sx2) - (Sx * Sx)) 
End If 
End If 
End Function 

Sub ReaO 

'calculates and records Rea data 

' f ind index number (m) of the gal . 
'center, the center already chosen 
redom: 

IData(O) = 

For k = 1 To p 

If IData(O) < IData(k) And 
IData(k) <= dmax Then 



IData(O) = IData(k) 
AData(O) = AData(k) 
m = k 
End If 

Next 

For k = 1 To p 

If k <= m And AData(k) > Then 

AData(k) = -AData(k) 
End If 

If k > m And AData(k) < Then 

AData(k) = -AData(k) 
End If 

Next 

're-record data as surface brightness 
'and angular data in arc seconds 
Datal .Refresh 
Datal . Recordset .MoveFirst 
For k = 1 To p 

Datal . Recordset . Edit 
Datal . Recordset . Fields ( " INDEX" ) 
.Value = CSng(k) 

If IData(k) < 1 Then 
Datal 

. Recordset . Fields ( "VALUE" ) _ 
.Value = CSng(25) 
Else 

Datal 

. Recordset . Fields ( "VALUE" ) _ 
.Value = CSng(VegaZP - 2.5 * _ 
(Log(CDbl(IData(k)) / ExpTime / Area)_ 
/ Log(10#))) 
End If 
Datal 

. Recordset . Fields ( "DEGREES " ) _ 
.Value = CSng(36G0 * AData(k)) 
Datal .Recordset .Update 
Datal .Recordset .MoveNext 
Next 

Datal . Recordset . Close 

'pos going 
'find starting k 
k = 

Do While IData(m + k) > 0.95 _ 
* IData(m) 

k = k + 1 

Loop 

If IData(m + k) < 0.8 * IData(m_ 
) Then k = -1 

CountH = 

TSlope = -10000000 

Sb = -20000000 
redoSlopePS: 

SbG = Slope2(G) 'reset slope 
' numbers 

For j = 1 To 2 'initialize 
' slope 

k = k + 1 
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If IDataCm + k) < _ 
IDataCm + k - 1) Or k = Then 
SbO = Slope2(m + k) 

Else 

j = j - 1 
End If 

Next 

CountL = 

Do While CountL < 10 And_ 
IDataCm + k) > 50 And (m + k + 2) > 1 
k = k + 1 

If IDataCm + k) > 5 Then 
'for false data pixel 

Sb = Slope2Cm + k) 
If SbO < Sb Then 

CountL = CountL + 1 
CountH = 
SbO = Sb 

Else 

CountH = CountH + 1 
CountL = 
SbO = Sb 
End If 

Else 

IDataCm + k) = 60 
End If 

Loop 

Rip = m + k - CountL 

'neg going 
'find starting k 
k = 

Do While IDataCm + k) > 0.95 _ 
* IDataCm) 

k = k - 1 

Loop 

If IDataCm + k) < 0.8 * IDataCm) _ 
Then k = 1 

CountH = 

TSlope = 10000000 

Sb = 20000000 
redoSlopenS : 

SbO = Slope2C0) 'reset slope 
' numbers 

For j = 1 To 2 'initialize 
' slope 

k = k - 1 

If IDataCm + k) < _ 
IDataCm + k + 1) Or k = Then 
SbO = Slope2Cm + k) 

Else 

j = j - 1 
End If 

Next 

CountH = 

Do While CountH < 10 And 
IDataCm + k) > 50 And Cm + k - 2) > 1 
k = k - 1 



If IDataCm + k) > 5 Then 
'for false data pixel 

Sb = Slope2Cm + k) 
If SbO < Sb Then 

CountL = CountL + 1 
CountH = 
SbO = Sb 

Else 

CountH = CountH + 1 
CountL = 
SbO = Sb 
End If 

Else 

IDataCm + k) = 60 
End If 

Loop 

Rln = m + k + CountH 

Real = AbsCADataCRlp) - ADataCRln)) 
If AbsClDataCRlp +2)) < 10 Then 

Real = AbsCADataCm) -_ 
ADataCRln)) * 2 

IDataCRlp) = 
End If 

If AbsClDataCRln)) < 10 Then 
Real = AbsCADataCm) -_ 
ADataCRlp)) * 2 

IDataCRln) = 

End If 

'record found r_a in Data2 which 
'is Ein Excel file 

Data2 .Refresh 

Data2 . Recordset . MoveFirst 

Do While Data2 
. Recordset. Fields C "USED") .Value = "Y" 
Data2 . Recordset . MoveMext 

Loop 

If Data2 . Recordset . Fields C "USED" ) _ 
.Value <> "Y" Then 

Data2 . Recordset . Edit 
If IDataCm) < 1 Then 

Data2 

. Recordset. FieldsC"PEAK") .Value = 0_ 
Else 

Data2 . Recordset 
. Fields C "PEAK") .Value = CSngCVegaZP - 
2.5 * CLogCCDblC IDataCm)) /_ 
ExpTime / Area) / LogCl0#))) 
End If 

If IDataCRln) < 1 Then 

Data2 . Recordset . Fields_ 
C"R1SBDPN") .Value = 
Else 

Data2 . Recordset . Fields_ 
C"R1SBDPN") .Value = CSngCVegaZP - 2.5_ 
* CLogCCDblClDataCRln)) / 
ExpTime / Area) / LogClO*))) 
End If 
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If IData(Rlp) < 1 Then 

Data2 . Recordset . Field.s_ 
("RISBDPP") .Value = 
Else 

Data2 . Recordset . Fields_ 
("RISBDPP") .Value = CSngCVegaZP -_ 
2.5 * (Log(CDbl(IData(Rlp)) /_ 
ExpTime / Area) / Log(10#))) 
End If 

Data2 . Recordset_ 
.Fields ("USED") .Value = "Y" 

Data2 . Recordset . Fields_ 
("SOURCE") .Value = Left(Filel_ 
.filename, Len(Filel. filename) - 4) 

Data2 . Recordset . Fields_ 
("2REA1") .Value = Real 

Data2 . Recordset . Fields_ 
("2REA2") .Value = Rea2 

Data2 . Recordset . Fields_ 
("PEAKDEG") .Value = AData(m) 

Data2 . Recordset . Fields_ 
("M") .Value = m 

Data2 . Recordset . Fields_ 



("RINIVALUE") .Value = AData(Rln) 

Data2 . Recordset . Fields_ 
("RIPIVALUE") .Value = AData(Rlp) 

Data2 . Recordset . Fields_ 
("MOPIX") .Value = p 

Data2 . Recordset . Fields_ 
("PA") .Value = PA 

Data2 . Recordset . Fields_ 
("INSTRUMENT") .Value = Instrument 

Data2 . Recordset . Fields_ 
("FILTER") .Value = Filter 

Data2 . Recordset . Fields_ 
("VEGAZP") .Value = VegaZP 

Data2. Recordset .Fields_ 
("EXPTIME") .Value = ExpTime 

Data2 .Recordset .Update 

Else 

Stop 'got a problem 
'with the data file 

End If 

Data2 . Recordset . Close 
End Sub 
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